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C^ ' Abstract 

We review the recent developments in the theory of the one-dimensional tight-binding 
Schrodinger equation for a class of deterministic ergodic potentials. In the typical examples 
the potentials are generated by substitutional sequences, like the Fibonacci or the Thue-Morse 
^ [ sequence. We concentrate on rigorous results which will be explained rather than proved. The 

\Q . necessary mathematical background is provided in the text. 
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1 Introduction 

The onc-dinicnsional Schrodinger equation 

{-(f/dx^ + V{x))i:{x) = E^P{x) 

has long been serving as a laboratory for studying a large variety of problems in Quantum Mechan- 
ics. Most often it appears as the eigenvalue equation of a one-electron energy operator, and then 
V{x) means the potential energy of the electron. In Atomic or Molecular Physics V{x) is typically 
fast-decaying or vanishes outside a bounded set, while in applications to Solid State Physics it is 
an extended function. Our concern will be the Schrodinger equation of Solid State Physics in the 
tight-binding approximation: the real axis is replaced by Z, the lattice of integers, and the differ- 
ential operator by the second order difference operator. The tight-binding Schrodinger equation 
describes an electron in a lattice in which the attraction of the ion cores slows down the interatomic 
motion of the electron: the kinetic energy has an a priori upper bound which does not exist in the 
continuous case. 

The potentials we are going to consider are on the midway between periodic and random; in 
physics literature they are sometimes called 'aperiodic' or 'disordered deterministic'. They all fall 
into the family of the so-called strictly ergodic sequences. A typical subfamily which appeared 
regularly in the lectures of this Winter School are the substitutional sequences as, for example, the 
Fibonacci sequence. As we will see later, the periodic and the uniformly almost periodic sequences 
are also strictly ergodic, but a random potential is usually 'only' ergodic. 

In this field physicists and mathematicians have been working parallelly during about the last 
fifteen years. The ambition of this course is to review and explain the contribution of the latter 
for an audience unfamiliar with the underlying mathematics. Most of the time, mathematicians 
consider the Schrodinger equation directly on the infinite chain Z or on the semi-infinite chain N. 
The equation is then the 'eigenvalue equation' of an infinite matrix which is named after Jacobi. 
If the diagonal (the potential) is an ergodic sequence, the spectrum happens not to be a discrete 
point spectrum. It is rather often - at least within the family of potentials we are going to consider 
- singularly continuous and is, as an ensemble of points, a (generalized) Cantor set. These notions 
and many others, necessary for the understanding of the specific results will be introduced in the 
forthcoming sections. 

We start, in Section 2, by the presentation of the main examples and the definition of strict 
ergodicity. The subject is developed in Sections 3, 4, 6 and 7. In Section 3 we give the basic 
ingredients of the theory of discrete Schrodinger equation in one dimension and introduce the 
notions of the transfer matrix and Lyapunov exponent. A subsection is devoted to the scattering 
problem and Landauer resistance. Section 4 contains some of the fundamental results for periodic 
potentials. Section 6 summarizes the general results obtained for one-dimensional Schrodinger 
operators with strictly ergodic potentials. Section 7 is a review of our actual knowledge about 
models with specific potentials, given by Sturmian and substitutional sequences. 

Section 5 is an introduction into spectral theory. It is far more the longest part of the course, 
with the definitions of spectrum, essential spectrum, eigenvalue, point and continuous spectrum, 
spectral projections, measures and spectral measures. Cantor sets etc., illustrated with examples. 
Subordinacy and the classification of the spectrum according to the asymptotic behavior of the 
solutions of the Schrodinger equation is also discussed here. 

The widely studied Almost-Mathieu or Harper equation is not presented in full detail but 
appears in examples throughout the whole text, and papers on it are included among the references. 

Some more or less simple facts are announced in the form of problems. The reader may only 
retain the statements or, if he wishes, check his understanding by proving them. The solutions can 
be found in Section 8. 

These notes are not written in a purely mathematical style. Proofs, unless they are really 
elementary, are replaced by explanations. This suffices for the physicist who wishes to understand 
the results but has no ambition to provide mathematical proofs of his own. The interested reader 
can find excellent monographs which present the subject with full mathematical rigor. Also, one 
should remember that nothing can replace the lecture of the original research papers. A rather 
detailed, although certainly incomplete, list of references is provided at the end of the notes. I 
thank Frangois Monti for helping me to construct this list. 



It is a pleasure to thank Frangois Dclyoii, Frangois Monti and Charles Pfistcr, who accepted to 
read the manuscript; their remarks brought considerable improvement to the final text. 

2 Main examples 

The tight-binding Schrodinger equation in one dimension is written as 

V'n-l + IJM+I + Vnlpn = Elpn (2.1) 

Here n runs over the integers, V — {Vn} is a real sequence (the potential), £■ is a real number (the 
energy), although sometimes it is useful to take it complex. The problem is to find the 'physically 
allowed' values of E, that is, those admitting at least one 'not too fast increasing' solution, and to 
describe these solutions. The following potentials will be considered: 

(1) Vn = V{uj)n — Acos27r(na + w), A ^ 0. The corresponding discrete Schrodinger equation 
is called the Almost-Mathieu equation (because the continuous version is the Mathieu equation). 
The case A = 2 is the Harper equation. 

(2) Vn = V{uj)n — \XA{na + ui), where \ ^ 0, A is the union of semi-open intervals of [0, 1) 
and Xa{x) = 1 if the fractional part of x falls into A and is zero otherwise. These are the so-called 
circle-potentials. An important special case is when A is a semi-open interval of length a: The 
{0, l}-valued sequences 

XiQ^^\(na + Lu) = [na + wj — [(n — l)a + lu\ 

X(a,a]{fT'a + Lu) = \na + lu'] — \{n — l)a + lu~\ (2-2) 

for irrational a are called Sturmian sequences. The notations [-J and [•] mean rounding downwards 
and upwards, respectively. If a = (yS — l)/2, we obtain the Fibonacci sequence. Generalized 
Fibonacci sequences (obtained through the substitution ^(a) == a"b and £,{b) = a with n > 1) are 
also Sturmian. According to an equivalent definition, Sturmian sequences are the sequences with 
minimal complexity among the non- ultimately-periodic sequences (see e.g. [1], [23]). 

(3) F is a substitutional sequence: 

Let ^ be a primitive substitution ([114]) on a finite alphabet A = {oi, ...,0^}, / a non-constant 
real function on A. Suppose that ^(a) = ...a and ^(b) — b... for some a,b Cz A and let uj = uv where 
u and V arc the left- and right-sided fixed points of ^: u — ^°°(a) = ...a and v = C°°(&) = b.... 
Then V{uj)n = fi'-^n)- A comment on the concatenation uv will follow in Section 6. Generalized 
Fibonacci sequences have a natural definition on Z via the formula for Sturmian sequences. 

The potentials defined in (l)-(3) are strictly ergodic, that is, minimal and uniquely ergodic: 

(i) V is minimal. 
Let T be the left shift, and denote V{Tuj) the shifted potential: V{Tu!)n — V{uj)n+i- The notation 
(T applied to lo) may seem curious, however, in all cases Tlo can be given a well-defined meaning. 
For potentials (1) and (2) 

Tcu = LO + a (mod 1) 

In case (3), Tuj is the left-shifted sequence, 

{Tuj)n = UJn+1 

The orbit of V{ljj), i.e., the sequence of translates, V{T'^(jj) for k = ..., —1,0, 1, ..., has pointwise 
convergent subsequences: W = {Wn} is called a pointwise limit of translates of V{ui) if for a 
suitably chosen sequence {h}, V{T'^'uj)n tends to Wn for every n. Clearly, W — V{lj'), where uj' 
is the limit of T'^^uj in cases (1) and (2) and is a pointwise limit of translates of uj in case (3). 

The orbit of V together with all its pointwise limits is called the hull of V and is denoted 
by ^(V). This is a kind of closure of the set of translates. The hull can be constructed for any 
infinite sequence. If V is periodic with period length L, the hull is a set of L elements. The 
potentials of (1) and (2) are periodic if a is rational. If a is irrational, the hull oi V = V(oj — 0) 
is {y(w)}(jg[o,i)! This is the typical situation: the hull is much larger (uncountable) than the set 
of translates (countable). If 5* is an arbitrary sequence and S" is a pointwise limit of translates of 



S, the hull of S" is usually different from the hull of S. (Example: Sn ^ I — 1/n, n = 1, 2, .... The 
translates have a single pointwise limit, S'^ = 1. The only element of 17(5") is S" itself.) 

The potentials defined in (l)-(3) have the property that for any W G 0(F), V,{W) = ^{V). 
This is what the minimality of V means. The V^-independent hull is denoted by fl. Clearly, every 
element of fl is minimal and therefore the 'flow' (J7, T) itself can be called minimal: every trajectory 
is dense in 17 (its closure is il) or equivalently, the only closed and shift-invariant subsets of fl are 
the empty set and fl. 

Problem 1. (Cf. Queffelec [113]) Let s ~ {sn}-oo be a bounded complex sequence, VL{s) its 
huU, 

where the bar means closure with respect to pointwise convergence. The following are equivalent. 

1. s is minimal, i.e., for every t € f7(s), f7(t) ~ il(s). 

2. For every t e ri(s), s e VL(t). 

3. s is almost-periodic in the metric 

oo 
n— — oo 

i.e., for each e > there exists an f^ < oo and a sequence {n/^}^^^^ of integers with gaps bounded 
by ^g, < nk-\-i — Tik < ^ei such that for all k 

oo 
n— — oo 

A sequence is called recurrent if it is the pointwise limit of its own translates. A minimal 
sequence is recurrent but a recurrent sequence may not be minimal. A counterexample is provided 
by the sequence of digits 1234567891011121314..., obtained by concatenating the positive integers, 
written, e.g., in base 10. It is the pointwise limit of its own left shifts, but is not minimal: words 
are repeated with increasing gaps. The hull of this sequence is the set of all sequences composed 
of the digits 0, 1, ...,9. 

The different meanings, in different cases, of the label uj should not disturb the reader. In (1) 
and (2), w is a number of the interval [0, 1); in (3) it is an element of the hull of uv, and uv is also 
a minimal sequence. 

(ii) V is uniquely ergodic. 
Recall that the dynamical system (il, T, p), where p is a T-invariant probability measure, is ergodic 
if for any p-integrable function / 

N 

for p-almost every uq (i.e., apart from a set of ujo's of zero p-mcasure). Choosing / to be the 
characteristic function of a set A, it is seen that ergodicity implies that (in fact, is equivalent to) 
p{A) = or 1 for every T-invariant set A. 

The unique ergodicity of V means that there exists a unique shift-invariant probability measure 
p on n{V) {p{T^'^A) = p{A) for any measurable A C il). 

Problem 2. If s is a uniquely ergodic sequence with shift-invariant probability measure p then 
(fl(s),T,p) is ergodic. 

In case (1) (Almost-Mathieu) , the most natural way to obtain this measure is to define, for 
<a< 6< 1, 

p{Aab)=b-a , whCTC Aab ^ {V{uj)}a<uj<b (2.3) 

and to extend p to countable unions and intersections of the sets Aab by using the properties of 
measures. 

For substitutional potentials the usual construction is to assign probabilities to the so-called 
cylinder sets and to extend p to countable unions and intersections of cylinders, liv — {vi,V2, ■■■,Vk) 



is a finite collection of real numbers and n is an integer, the cylinder [v, n] is a family of potentials 
taken from the hull, 

[v, n]^{V en■.Vn+^^v„^ = l,..., k} (2.4) 

Noticing that the relative frequency of the 'word' v is the same in all V E H., the probability 
p{[v, n]) can be defined as this relative frequency. 

For circle potentials a shift-invariant p can be derived naturally either through Eq. (2.3) or 
through word frequencies. For the Almost-Mathieu potential cylinders can also be defined (by 
replacing, in Eq. (2.4), Vn+i = Vi by Vn+i € li, where li is an interval), and one can obtain a 
T-invariant probability measure through them. Uniqueness implies that the different constructions 
yield the same measurable sets and measures. 

Physicists never worry (do not even have to know) about the hull; mathematicians all the time 
struggle with it. The reason is that in laboratory or computer experiments one deals with finite 
samples, and these may come from any element of the hull: therefore it is sufficient to keep just 
one. (However, the physicist's sample averaging is nothing else than averaging over the hull with 
the measure p.) The mathematician is naturally led to the notions of hull and strict ergodicity 
by working on the infinite system. He (she) has to consider 1/ as a p-distributcd random variable 
on a.. The results she (he) obtains are valid either for a given V G ft oi for every V £ ft oi, 
most often, for p-a.e. (almost every) V € ft, i.e., p{{V G ft: the result is not vahd}) — 0. Notice 
that the distinction between 'every' and 'almost every' is by no means harmless. Indeed, because 
p is a continuous probability measure (that is, p({V^}) = for every one-point subset of fl), any 
countably infinite subset of the hull has zero probability. A surprising new result about an even 
larger (uncountable) set of vanishing p-measure, for which a claim (localization) is not valid, will 
be quoted in Section 5.2. 

Problem 3. If p is a shift-invariant probability measure on ft and the trajectory of every point 
is infinite then p is continuous. 

3 General results on the Schrodinger difference equation 

3.1 Basic observations 

The first two terms of Eq.(2.1) correspond to the action of the kinetic energy operator on the wave 
function, and any physicist would prefer to replace them by 

This is, however, a difference of no importance: Omitting 2i/'„ means shifting the zero of the energy 
by 2, and changing the sign of the first two entries in Eq.(2.1) amounts to a unitary transformation, 
thus leaving the spectrum of the energy operator unchanged. Therefore a global change of sign 
of V is also of no importance: we get the mirror image of the spectrum. This simple fact can be 
expressed in a more shocking form: If V gives rise to a localized state at energy E, —V will give 
rise to a localized state at energy ~E. Or still in another form: There can exist extended states 
at energies everywhere below the potential, E < Vn for all n. 

Without saying otherwise, we always consider the Schrodinger equation on the infinite lattice 
Z. Equation (2.1) is a homogeneous linear second order difference equation. As a consequence, for 
any (complex) E the solutions form a two-dimensional linear vector space: 

• If Tp^ and f/;^ solve (2.1), then xp — ciip^ + C2ip'^ solves it for any complex ci and C2. 

• There exist two linearly independent solutions tp'^ and ip'^. 

Two nontrivial (^ 0) solutions are linearly dependent if and only if they are constant multiples of 
each other. When saying that Eq.(2.1) has a unique solution of some property, we will understand 
'unique linearly independent'. 

Writing (2.1) for Tp^ and V'^ (not necessarily linearly independent), multiplying the first equation 
by f/'n and the second one by 'ipn a-^d subtracting we find that 

■(Pl+i^l - i^l^l+i = V'iV'Li - V-i-iV-^ = • • • = ipl^l - ipl^l 

This constant is the Wronskian and will be denoted by W[ip^,^p'^]. As a consequence, 



• ip'^ and ip"^ are linearly independent solutions if and only if Wlip^jip"^] 7^ 0. 

• There can be at most one solution going to zero as n — > 00. If i/i^ is such a solution, then for 
any other (linearly independent) solution ip 

\^nf + l^ri+iP ^^ 00 as n -> 00 

The same conclusion holds for the limit n — > —00. As a consequence, there can be at most 
one solution decaying on both sides, and if such a solution exists, all the others blow up on 
both sides. 

3.2 Transfer matrices 

We may rewrite the Schrodinger equation in the redundant vectorial form 



where 

By iterating, 
and similarly. 



*" - ( '^;t 1 (^-1) 



*„ = r„r„_i • • • Ti*o - Ti^„*o 



The matrices T are called transfer matrices. Their determinant is 1: this is clear for T„ and its 
inverse and holds for the others because the determinant of a product of matrices is the product 
of the determinants. 

The determinant of Ti^^ can be seen as the Wronskian of the two 'standard' solutions ijj^ and 
ip"^ with initial conditions 

V-i = 1 Q^ipl 
In terms of these solutions the transfer matrix can be written as 

showing that indeed, 

detri^„ = i^[V'\V''] = i 

The general form of the transfer matrix from site m to site n is 

T,n^n =A=( ^j , ad-bc=l (3.3) 

Here a, b, c, d are polynomials of E with real coefficients and of degree < |n -— to| + 1. In particular, 

trA = a + d = Sl"-™l+i + --- 

For real energies the transfer matrices are elements of the group SL(2,R). The characteristic 
equation of a matrix A in this group reads 

det(A - A) = A^ - A tr A + 1 = 

Therefore the two roots satisfy the relations A1A2 — 1 and Ai + A2 — ti A. We can distinguish four 
cases: 



1. The two roots are real, |Ai| = |A^"'^| > 1, the corresponding eigenvectors are also real and 
nonorthogonal, in general. |tr A| > 2, the matrix is hyperbolic. 
For example, the 'hyperbolic rotations' 

cosh X sinh x 
sinh X cosh x 

are special (symmetric) hyperbolic matrices. Ai_2 — exp(±x), the corresponding eigenvectors 
are the transposed of (1 ±1). 

2. Ai = Aj = e**^. If 7^ mr. There are two complex eigenvectors, tr A ~ 2cos(^, the matrix is 
elliptic. As an example, the ordinary rotation matrix 

cos X — sin X 
sin X cos X 

is elliptic. Ai,2 = exp(±za;), the eigenvectors are the transposed of (1 =F i). 

3. Ai := A2 = ±1. There is a unique, real, eigenvector (the geometric multiplicity is 1). |tryl| = 
2, the matrix is parabolic. An example is 

1 b 
1 

with b y^ 0. The unique eigenvector is the transposed of (1 0). 
A. A ~ ±1 where / is the unit matrix. 
Some assertions will concern the norm of transfer matrices. The natural matrix norm 

||yl|| = sup||Am||/||u|| 
associated with the vector norm ||u|| = (|uip + |u2p)^ satisfies 

1 o /I, 



\\Ar = max. eigenvalue of A* A = ^\\AWi + [ -\\A\\t - 1 ) = P'^f > 1 



(3.4) 



Here ||A||2 is the trace norm, 

\\A\\l = tTA*A = a^+b^ + c^+d'^ = \\A-^\\l >2 (3.5) 

(cf. Eq.3.3), which proves also the last inequality in (3.4). 

Let us return to the trace of the transfer matrix as a function of E. 

trri^„(£;)-s" + ---=p„(£;) 

is an nth degree polynomial with n distinct real roots and the property that 

^PniE) ^ if \pniE)\ < 2 

In Figure la we represented the generic situation for n = 5: There are n disjoint intervals of 
the energy within which the transfer matrix is elliptic, separated by intervals of hyperbolicity. 
Figure lb shows an atypical case for n ~ 4: Some intervals of hyperbolicity disappear because the 
derivative of the trace vanishes at the same E where the trace takes on the value 2 or —2. 

Problem 1. For A eSL(2,R), show that |try4''| < 2 if and only if |tr^| < 2. Plot trr3_^„(S) 
on the top of Figure 1. 



3.3 Lyapunov exponent 

Consider 



In = -p-rlnllTi^ 

\n\ 



This number is nonnegative (because the norm is > 1) and depends on E and the parameters of 
the potential. From 7„ we obtain, as n — > ±00, 

7± = Um sup 7„ , 7 ± = hm inf 7„ 

We find the same numbers by using the trace- or any other norm. In fortunate cases the hmits 
exist (i.e., limsup and liniinf coincide) and give 7± = ^±{E^V). These we call the Lyapunov 
exponents. If 7+ =7-, the common value is denoted by 7(i?, V). 
To elucidate the meaning of the Lyapunov exponent, recall that 

cf. equations (3.1), (3.2) and (3.5). Let, e.g., n — > +00 and suppose that 7-|_ exists. If 7-1- = 0, 
every solution is subexponential. Suppose that 7+ > 0. Then asymptotically 

In 11*^11 xn7+ 

for i = 1 or 2 or both and, since any solution is a linear combination of the two standard ones, the 
same asymptotics is valid for all solutions, except perhaps a unique one (which is asymptotically 
smaller). The Multiplicative Ergodic Theorem due to Osceledec [104] and Ruelle [123] asserts that 
there is indeed a unique solution ip with an exponential decay: 

ln||*„|| X -n7+ 

If both of 7± are positive, it may happen that the unique solution decaying exponentially to the 
left coincides with the unique solution decaying exponentially to the right: we have an exponentially 
localized state. If this is not the case, for the given energy all solutions increase exponentially at 
least in one direction. Such an energy is considered as 'physically forbidden' and the solution 
physically irrelevant, because there can be no electronic state inside the solid, corresponding to it. 

Let us make this statement somewhat more explicit. We may ask: Provided that the elec- 
tron is in the interval [—L^L], what is the probability of finding it at site n? The answer is 
IV'nj^/ X]fe=-L I'/'fej^- This number is of the order of unity close to L or —L and exponentially 
decaying towards the inside. So in any finite sample, at the given energy, the electron cannot 
penetrate into the bulk, it is exponentially localized at the boundaries. Also, one can learn from 
generalized eigenfunction expansion that solutions increasing faster than any power of |n| will not 
appear in the expansion of wave packets or computation of matrix elements and are therefore of 
no relevance to Physics. 

3.4 Scattering problem: Landauer resistance 

Suppose we have a finite sample: a potential given at n = 1, ..., L. Our aim is to discuss scattering 
of an electron of energy E on the sample. We devise our scattering experiment so as to have access 
to all energies. To this end, the potential is extended to Z with constant values, 

_ J cji if n < 
1 CJ2 if n > L 

chosen so that \uji — CJ2I < 4. Then the scattering problem is defined for energies satisfying the 
inequalities 

\E-ui\<2 , \E-U2\ <2 

(Notice that E < uji is allowed, see the remark in the first paragraph of Section 3.1.) The general 
solution has the form 



^n - { (j^k^n ^ £)e-*fe2«'^ ^ > ^ (3-6) 



(the plane waves extend to n = 1 and n = L) where the real wave numbers ki are related to E and 
LUi via the equations 

2 cos ki ~ E — bJi , i ~ 1,2 

There is no loss of generality, if we suppose that < fci < tt. With this convention, the scattering 
problem amounts to set A = 1, B = r, C = t and D = and to solve for t and r. The Landauer 
resistance of the sample of length L is defined as Rl = |rp/|tp. 

Tp with the scattering boundary condition is a complex solution of the Schrodinger equation. 
Therefore, its complex conjugate is a linearly independent solution to the same (real) energy. The 
Wronskian between tp and ip* is a constant times the current density: its constancy physically 
means current conservation. Compute the Wronskian at n = and ed, n ^ L. One finds 

W[ij,^*] =2i(l-|r|2)sin/ci =2i|tpsinfc2 

so that 

|<|^sinfc2 = (1 - |r|^)sinfci 

This relation is well-known in the special case of fci — k2. Because the sines are positive, \r\ cannot 
exceed 1. 

We want to express Rl in terms of the elements a, 6, c, d of the transfer matrix Ti^^. (Here 
we use the notations of Eq.(3.3)). This can be done directly by solving the linear problem 

ie*'=2i J ^ y c d ) \ 1 

but there is a more elegant way to do it. We can express the Landauer resistance also in terms 
of another type of transfer matrix, usually applied for the Schrodinger equation in the continuum. 
This transfer matrix, 

' a 13 ' 



7 
connects the coefficients A and i? to C and D, sec Eq.(3.6): 

In contrast with Ti^l, Y depends on ki and k2. Let now ip denote the solution with A — 1 and 
-B = 0; it follows that C — a and D = 7. tp*, the complex conjugate of tp, is also a solution with 
^ = and B — 1 and, thus, with C = /3 and D ~ 5. We find therefore 

7 = /?* , 5 — a* 

Furthermore, the Wronskian of tp and tp* computed at n = and aX n — L yields 

W[tP,tP*] =2isin/ci =2i(|ap- |/3|2)sinA:2 
and thus 



The scattering problem 



trivially gives 



from which 



dety = \aY - \li\ = sinki/ s\n.k2 



t \ y/ 1 



/ \ r 



r^-fi a , t^—^—- (3.7) 

a sm K2 



Rl — 



2sin^fc2 
sin fci 



10 



Another expression of Rl in terms of the trace norm of Y can also be obtained. 

sinfc2 

With this we find 

1 ,, ,,oSin fco 1 sinfco 

'^"4ll'^ll2^in2fci 2sinfci 

These formulas are frequently used with fci — k2, see e.g. [46], [85], [82]. It remains to connect the 
two different kinds of transfer matrices. This can be done by simple inspection and results 

Solving this equation for a and /3 and substituting into Eq.(3.7) we obtain 

t = -2i sin/cie-*'=2i(ag-ifci _^ 5 _ e^'=2(gg-ifci ^ ^))-i (3 g) 

r^- "^, ^: ^, +^^^, (3.9) 

and 

RL{ki,k2) = (4sin2fci)-i|ae^'=i + 6 - ce*('=i+'=2) - de^'^^p (3.10) 

Recall that R^ depends on E via a, 6, c, d which are real polynomials of E. The ordinary choice 
for the potential outside the sample is wi = u;2 = 0. With this we get fci = A:2 = fc, 2cosA: ~ E 
and can 'test the sample' at energies between —2 and 2. 

A nice and simple form is obtained for ki = k^ = 7r/2 which corresponds to wi = W2 = E: 

Rl{^/2, V2) = |za + 6 + c - id\''/A = (UTi^Llla - 2)/4 

A truly intrinsic resistance is the minimum of Rl over all fci , k2 which, in general, cannot be 
given in a closed form. The minimum over fc2, when fci is fixed at 7r/2, is still easy to find: 

mini?L (71/2,^2) - i?L(V2,arg^^) - {l/A)Wa^ + h^-^c^ + cPf 
fc2 c — id 

If the Lyapunov exponent is positive, the asymptotic behavior of the resistance is the same for all 
fci, ^2: Rl grows with L like exp{27+i}. 

4 Schrodinger equation with periodic potentials 

Periodic differential equations have a fully developed beautiful theory (Floquet theory, see e.g. 
[45]) which can be adapted without any difficulty to the case of difference equations. Moreover, 
strictly ergodic sequences are well approximated by periodic ones, and most of the results discussed 
in Sections 6 and 7 deeply rely upon this fact. In this section we therefore survey the periodic case. 
Let the potential be periodic with period L > I finite: 

Vn+L = Vn for any integer n 

Consequently, 

J-rn+L^n+L ^J-rn^n j S-h 171,11 

We exploit this by writing n = kL + m where k — k{n) and < m = m{n) < L arc uniquely 
determined and, e.g. for n > 0, 
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According to the value of E we can distinguish several cases (see also Figure 1): 
(1) The transfer matrix is elliptic or hyperbolic, 

|trTi^i(S)|-|A + A-i|^2 

where A and A^^ are the eigenvalues. In this case there exists a matrix S (generally nonunitary) 

such that 

A 
A- 



S Ti^lS — I r, , _i 



Therefore wc can write 



The initial condition 



yields the solution 



while with 



*„ = ri^„5( -; \ ) 5-1*0 



A*^ 
A- 



*o = 5( J 



*n = ^ Ti^jn(n)S 



*0 = 5> ^ 







we get 



^2 ^ A-'=(")Ti^„(„)5 (^ J ) 



Because m{n) is an L-periodic function, A**^ in $i'^ multiplies an L-periodic function of n. There 
are two cases: 

(1.1) The transfer matrix is elliptic, 

\ — e^^ , tr Ti^ l{E) = 2 cos rj , r/ 7^ integer x tt 

The two solutions cf)^''^ are Bloch waves and E is in one of the (at most L) 'stability intervals'. 
The Lyapunov exponent exists, 7+ = 7- = "f{E) = 0. The Landauer resistance _R„ is a bounded 
quasiperiodic function of n (periodic, if ry is a rational multiple of tt). 

(1.2) The transfer matrix is hyperbolic, say, |A| > 1. The two solutions can be written as 

0i=AL"/^Jp„ , 02^A-L"AJ<z„ 

with L-periodic sequences Pn,qn- The energy is in an open - maybe semi- infinite - interval. The 
finite intervals are called 'forbidden gaps', they separate the stability intervals. The Lyapunov 
exponent "f{E) =■ L~^ In |A|. The Landauer resistance i?„ grows exponentially (with exponent 27). 
(2) The transfer matrix is parabolic. There exists a unitary matrix U such that 

The solution can therefore be written as 

/ 1 kh \ 
*„ = Ti^^U (01) '^''^o 

The initial condition 



gives rise to the periodic solution 
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and any other solution oscillates and grows linearly with a slope h/ L. The energy is a bound- 
ary point of a stability interval. The Lyapunov exponent is zero, the Landauer resistance has a 
quadratic increase. 
(3) Accidentally it may happen that Ti^l(E) = ±7. Then the general solution 

is periodic for all initial conditions. This situation is quoted as the 'closure of a forbidden gap' (cf. 
Figure lb). In most concrete cases, it is difficult to say whether it occurs or not. If all the gaps 
are open, there are L — 1 of them, separating L disjoint stability intervals. Some examples : 

(i) For the Almost-Mathieu equation with a — p/q rational and qLo/iT noninteger all the (q — 1) 
gaps are open. If cj = 0, p is odd and q ~ 4to, then one gap is closed (Bcllissard, Simon [21], Choi, 
Elliot, Yui [31]). 

(ii) For circle potentials (2.2) with rational a all the gaps are open. 

(iii) We can obtain periodic potentials by iterating a substitution on a letter a finite number 
of times and repeating the result periodically. Doing this with the period-doubling substitution 
{(T{a) = ab, a(b) = aa), for the resulting equation all the gaps are open (Bellissard, Bovier, Ghez 
[15], [27]). For the Thue-Morse substitution (cr(a) ~ ah, a{b) — ba) the trace map ([9]) 



t'n+2 



- 2 = (.T„+i - 2)a;2 (n > 1) 



with Xn — tr Ti^2" reveals that a;„ is a 2" degree polynomial of E, and a;„+2 — 2 has double zeros 
at the 2" simple zeros of Xn- This shows (see Fig. 1) that at least one fourth of the gaps is closed! 
In fact, 1/4 is the precise proportion of the closed gaps ([12]; see also [10] for the corresponding 
phonon problem). 

(iv) For the best periodic approximations of a hierarchical potential all the gaps are open (Kunz, 
Livi, Sfito [90]). 

We close this section with the conclusion that the physically relevant solutions belong to energies 
which fall into the stability intervals: 

\tTT,^L{E)\ < 2 

These solutions are the well-known Bloch waves. The stability intervals are nothing else than the 
bands in the spectrum of the periodic Schrodinger operator. 

5 Spectral theory 

Let us come back to the question we asked already in the previous sections. The Schrodinger 

equation on R"^ or Z'^ can be solved for any E. Which are the physically relevant solutions? What 

can be said about the corresponding energies? Spectral theory provides the following answer: 

(i) To compute matrix elements of functions of the energy operator, it suffices to consider only the 

polynomially bounded solutions of the Schrodinger equation. 

(ii) Every subexponential solution (growing slower than any exponential) belongs to an energy 

which is in the spectrum of the energy operator. 

5.1 Schrodinger operator, £^-space and spectrum 

Define an operator H on the sequences ip — {ipn} by 

Then Eq.(2.1) is equivalent to Hip ~ Etp. We call H the Schrodinger operator; it can be represented 
with the infinite tridiagonal (Jacobi) matrix (denoted by the same letter) 



H 



/ . 






\ 




1 K-1 


1 






1 


Vn 1 




V 






• / 



(5.1) 
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Consider 

l^ n— — oo ) 

This is a separable Hilbert space (it contains a countable dense set) with the inner product {(p, tjj) = 
S^-oo 'fni'rn and i7 is a selfadjoint operator on it. In what follows, we use the notation T-L for a 
separable Hilbert space. We call the elements of T-L vectors. The norm of a vector -0 is defined by 
II V'lP = (V'l V')- It should not be confounded with the norm of the C^ vectors which we introduced 
in Eq.(3.1) and denote by capital Greek letters. 

Suppose, for the sake of simplicity, that F is a bounded sequence. The spectrum a{H) of H 
is the set of E values for which E — H has no bounded inverse: One can find some vector ^ such 
that the equation 

{E-H)ip = xP (5.2) 

has no vector solution for ip (but may have a solution outside "H, see Problem 5.2.2) or it has more 
than one vector solutions. In the second case (E — Hjip = has a nontrivial vector solution, that 
we call an eigenvector of H, belonging to the eigenvalue E. (Recall that the spectrum of a finite 
matrix A is the set of its eigenvalues, and 

det(A — A) = if and only if (A — A)^^ does not exist.) 

The spectrum is real, nonempty and closed, i.e., contains all of its limit points. 
Example. Let us decompose H into off-diagonal and diagonal parts, H — Hq + V. Then 
a{Ho) = [-2,2]. If V is L-periodic, 

a{H)^{E:\tTT,^L{E)\<2} 

In either cases, the Schrodinger equation has no vector solution, thus there is no eigenvalue in the 
spectrum. 

If E is an isolated point in cr(-ff), then E is an eigenvalue. 

The essential spectrum, (JcssiH), of H is the set of all non-isolated points of (y{H), together 
with the isolated eigenvalues of infinite multiplicity (do not occur in one dimension) . Clearly, this 
is a closed set. An important characterization of the essential spectrum is due to Wcyl: 
E is in crcss(^) if and only if it is an approximate eigenvalue of H in the following sense. One can 
find a sequence (/?" of vectors such that ||(/5"|| — 1, </?" goes to zero weakly (in ^^(Z) this means 
that (/jJJ — > for each k fixed: either (/j" spreads over larger and larger domains or it remains well 
localized but 'escapes' to infinity) and 

H<^" = S</j" + e" (5.3) 

where e" are vectors, ||e"|j — )■ 0. 

If is not an eigenvector but it is a subexponential solution of Hijj = Eip, one can construct 
such a 'Weyl sequence', so E Cz crcss(-ff)- If ip is polynomially bounded, we call it a generalized 
eigenvector (eigenfunction) of H, and E a generalized eigenvalue. 

Problem 1. Construct a Wcyl sequence for a subexponential solution tp which is not an eigen- 
vector. Hint: Replace ipk by for fc < and k > n and normalize. 

5.2 Point spectrum 

The point spectrum of H is the set of all eigenvalues of H. It is denoted by app{H). The point 
spectrum may be empty, finite or countably infinite. The pure point subspace of H is the subspace 
spanned by all the eigenvectors of H. It is denoted by "Hpp. 

We say that H has a pure point spectrum, if its eigenvectors form a basis in the Hilbert space, 
that is, H = "Hpp. In this case the spectrum is the closure of the set of eigenvalues: 



a{H) = app{H) = app{H) U {hmit points of o-pp(_ff)} 

Example. A nontrivial counterexample known for every physicist is the hydrogen atom. Here 
the Hilbert space is L^(R^). The energy operator of the hydrogen atom has infinitely many bound 
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states (eigenvectors) which, however, do not span %. So "Hpp is an actual part of %. CTpp is a 
discrete set of negative numbers with a unique accumulation point, 0. Any E > Q belongs to the 
spectrum but is not an eigenvalue, only a generalized eigenvalue. The corresponding generalized 
eigenvectors are bounded: they are products of a plane wave with a confluent hypergeometric 
function. The generalized eigenvectors ('scattering waves') can be used to build up wave packets, 
i.e., square integrable functions (vectors of %) which are orthogonal to every eigenvector. 

Problem, 1. Find a possible Weyl sequence to show explicitly that is in the essential spectrum 
of the energy operator of the hydrogen atom. 

(Tpp is called a dense point spectrum, if it is nonempty and has no isolated point. If, for a 
bounded recurrent potential, there is a point spectrum, it is a dense point spectrum (sec Problem 
6.2.2). 
Examples of dense pure point spectra 

1. Let H = V, i.e., a diagonal matrix with Vn ~ cos27rriQ;, where a is irrational. V has a pure 
point spectrum, dense in [—1, 1]: 

app(y) = {F„}^.o , <j{v) = [-i,i] 

The eigenvector belonging to the eigenvalue Vn is (5", the unit vector concentrated on the site n. 
The set of eigenvectors is the canonical basis in i?^(Z). 

Problem 2. Choose E in [—1,1] but E yi^ Vn, any n. Use the definition of the spectrum to show 
that E G a{V). 

2. Anderson localization. Let Vn be identically distributed independent random variables. Let 
V{uj) = {V{u!)n} denote a realization, H{lu) = Hq + V{ijj). 

(i) Pastur [106]: There exists a set E cR such that 

• a{H{uj)) ~ E with probability 1. 

• E contains no isolated points. 

• Any _E G E is not an eigenvalue with probability 1. 

(ii) Kunz, Souillard [91] (simplified): Let the probability distribution r{x) of the T^i's be con- 
tinuous, and r{x) 7^ if and only if a < a; < 6. Then 

• the spectrum of H{uj) is pure point with probability 1. 

• E = [a - 2, 6 + 2] 

• All the eigenvectors are exponentially localized. 

(iii) Carmona, Klein, Martinelli [28]: Let Vn be Bernoulli-distributed, 

with probability p 



" ^ b with probability 1 — p 
Then the spectrum is pure point with probability 1, 

E = [-2,2]U[&-2,6 + 2] 

and all the eigenvectors are exponentially localized. 
3. Almost Mathieu equation (A > 0). 

Choose uj at random in the interval [0,1) (according to the uniform distribution). There is 
an appealing nonrigorous argument (Aubry- Andre duality, [4]), according to which the spectrum 
should be pure point for A > 2. This is not quite true, but seems to hold for good Diophantine a 
and almost every w: 

(i) If a is irrational, the spectrum is independent of u, a{H{uj)) = E(a, A). 

(ii) Sinai [131], Frohlich, Spencer and Wittwer [50]: If a is a good Diophantine number, namely, 
there exists some constant c > such that 

min{na — [na\ , \na\ — no} > cjr? for all n 7^ (5-4) 
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and A 3> 1, the spectrum of H(lu) is pure point for almost every w and the eigenvectors decay 

exponentiahy. Sinai's proof yields also that S(a, A) is a Cantor set (see below). 

Aubry, Andre [4], Thouless [138], [139]: The Lebesgue measure m of the spectrum is positive, 

, , ,, , , mCEia, A)) A — 2 

TO(E(a,A))>2|A-2|, I ^ '" > 



Eu — El A + 2 

{Ey: and Ei are the upper and lower boundary of the spectrum, respectively.) Recently it was 
shown by Last ([93], [94]) that for good Diophantine a, as for example in (5.4), 

m(S(a,A)) = 2|A-2| 

(iii) Jitomirskaya [69], [70] brought sensible improvements to the above results. She simplified 
the proof, weakened the Diophantine condition on a, put an explicit bound on A (A > 15), and 
showed that for all A > 2 and for almost every lo the closure of the point spectrum has the same 
Lebesgue measure as the spectrum itself. 

(iv) One may think that there is only some technical difficulty to extend the above results on 
localization from a.e w to every uj. This is not true: Recently, Jitomirskaya and Simon [71] proved 
that for any A > 2 and irrational a there is an uncountable set of w's which is dense in the interval 
[0, 1] (i.e., any point in [0, 1] is a limit point of this set) and for which the spectrum of i?(w) is 
purely singular continuous (see later). Clearly, this set is of zero Lebesgue measure, if for the given 
A and a localization occurs for a.e. w. This result provides an example of a singular continuous 
measure the support of which is a thick (positive-measured) Cantor set. 

5.3 Cantor sets 

Take a closed interval Cq of the real line. Cut off a finite number of open intervals which have no 
common boundary points with each other and with Cq. What remains is a closed set Ci without 
isolated points: a union of a finite number of closed intervals. Repeat the same procedure with 
Ci to obtain C2, and so on: continue it indefinitely, by following the rule that no interval is left 
untouched. Let C denote the resulting set. It has the following properties: 

1. C is closed. Indeed, each C„ is closed, Co D Ci D •••, therefore C — HJ^qCk, and any 
intersection of closed sets is closed. 

2. C is nonempty. Clearly, the boundary points c„i of C„ are not removed during the construc- 
tion, therefore C contains them for all n. 

3. There can be no interval in C, according to the construction rule. 

4. There is no isolated point in C . In fact, it is clear from the construction, that the set of 
points Cni is dense in C . 

A set with the above four properties is called a Cantor set. The boundary points of C"^, the 
complement of the Cantor set C, form a countable dense set in C, but the Cantor set itself is 
uncountable. C^ is an infinite union of open intervals, the closure of which is the whole real line. 
What is the Lebesgue measure of C? Suppose that m{Cn+i) — XniTi(Cn), then 

00 
m(C) = m(Co) Y\^ Xn 

n=0 

This number is greater than zero if and only if 

00 
y lna;„ > —00 

n=0 

which is true if and only if 

00 

y^(l -Xn) < 00 
n=0 

For Cantor's 'middle thirds' set a;„ — 2/3 and thus the Lebesgue measure is zero. 
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5.4 Continuous spectrum 

Let "Hcont = ^ppj the subspace orthogonal to "Hpp. By definition, H lias no eigenvector in "Hcont- 
If the only element of "Hcont is the null vector, the spectrum is pure point. If "Hcont contains a 
nonzero vector, H has (also) a continuous spectrum (JcontiH). Indeed, the restriction of H to 
^cont, denoted by H\'Hcont^ is a selfadjoint operator whose spectrum is nonempty, and 

<ycont{H) = cr(ii"|'Hcont) 

As a matter of fact, the dimension of "Hcont is cither or infinite ("Hcont is an _ff-invariant subspace 
[HT-Lcont C "Hcont), and H has eigenvectors in any finite dimensional invariant subspace) and in 
the second case crcont(-ff) is a nonempty closed set without isolated points. Thus we have the 
decomposition 

•H - -Hpp ffi -Hcont , a(i7) = crpp(iJ)Ucrco„t(i?) 

but notice that the closure of cTpp (even (Tpp itself) may overlap with CTcont • 

Just as to the points of app{H) there correspond the _ff-eigensubspaces of "Hpp, to certain 
uncountable subsets of (Jcont{H) there correspond infinite-dimensional iJ- invariant subspaces of 
^cont- The connection is made by spectral projections. 

5.5 Spectral projections 

If H has a pure point spectrum, it can be decomposed as 

Ee<jpp{H) 

where P{{E}) is the orthogonal projection to the eigensubspace belonging to the eigenvalue E (for 
the one dimensional Schrodinger operator all these subspaces are one dimensional) . The orthogonal 
projections are selfadjoint and idempotent, 

Pi{E})^ ^ Pi{E}) ^ Pi{E}r 
The orthogonality of the eigensubspaces belonging to different eigenvalues can be expressed as 

Pi{E})Pi{E'})^6E^E'Pi{E}) 
and, if the spectrum is pure point, 

^^ P{{E}) = / = identity 



E 



Beo-pp(H) 

Physicists ' notation: 
If _E is a nondegenerate eigenvalue, tp the normalized eigenvector, then 

P{{E}) =. W{^\ 

If E is degenerate. 



P{{E})^J2\^^)(^' 



where one has to sum over an orthonormal basis in the eigensubspace. 

In general, there is a fundamental relation between certain subsets of the spectrum and orthog- 
onal projections to iJ-invariant subspaces. Let A = (a, b] be a half-open interval of the real line, 
Xa its characteristic function, i.e., xa{E) = 1 if -E G A and otherwise. This function satisfies 

Xa{E) ^ XAiE)* = XA{Ef (5.5) 

One can define an operator xa (H) (formally obtained by substituting E with H) as a strong limit 
(see below) of polynomials oi H, which inherits the properties (5.5) of the real- valued function xa 
and is, therefore, an orthogonal projection; we denote it by P(A). Similarly, orthogonal projections 
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can be defined for more complicated sets, called Borel sets, as, for example, countable unions and 
intersections of half-open intervals (closed and open sets are Borel sets), and the family of the 
corresponding projections have nice algebraic properties: 

P(0) = , P(R) = / 

P(AinA2) = P(Ai)P(A2) 
and for Ai and A2 disjoint sets 

P(AiUA2) = P(Ai) + P(A2) 

The last equality holds also in a stronger form: If A^, i = 1, 2, ... is an infinite sequence of pairwise 
disjoint Borel sets then 

N 



Pi(JZi\)=s-\imJ2PiA, 



The sum converges in the strong sense, that is, for any ip ^Ti. 

N 

^P(A,)(^^P(U,^iA,)^ as A^^oo 

The above properties of P are characteristic to probability measures; the only difference is that 
P(A) are operators, so P is a projection- valued probability measure. This measure is called the 
(spectral) resolution of the identity, its values on Borel sets are the spectral projections. 
The spectral projections live on the spectrum of H: 

PiA) =P{Ana{H)) 

If P(A) 7^ 0, the linear subspace P{A)'H is nontrivial and invariant under H (because H commutes 
with P(A)). 

Consider now the converse relation: Let A be an _ff-invariant (closed) subspace and denote [A\ 
the orthogonal projection onto A. [A\ may not be a spectral projection. For example, if ip is an 
eigenvector belonging to a degenerate eigenvalue E and A is the one-dimensional subspace spanned 
by V', [A\ = |'0)(V'l is not a spectral projection. The smallest spectral projection larger than [A\ is 
P({-E}) (for projections Qi, (32, Qi < Q2 means QyH C Q2^)- Similar situation may occur in the 
continuous spectrum. The operator H is called multiplicity free if the orthogonal projection onto 
any iJ-invariant subspace is a spectral projection. Clearly, ["Hpp] and ["Hcont] are always spectral 
projections: 

[Hpp]=P(app(i7)) and [Hcont] - P(fT(iJ) \ ffpp (if)) 

Further important examples will be given later. 

For any vector ij) one can define a real positive measure /x^ by setting 

A.^(A) = (^,P(A)^) = ||P(A)^|p 

for Borel sets A. This is called the spectral measure associated with ■0. Spectral measures play an 
important role in the computation of averages and transition amplitudes: The matrix elements of 
functions of the energy operator can be obtained as integrals with respect to these measures. The 
spectral projection 

Pe = P{{-^.E]) 

is a monotonically increasing function of E in the sense that for Ei < E2 

Pe, - Pe, ^ P{{Ei, E2]) > 

(P > means {ip,Pip) > for all ip E H, which holds because {ip^Pip) ~ (Pip, Pip).) Therefore, 
with 

dPE = Pe+^e -Pe = P{{E, E + dE]) 
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one can write down the spectral decomposition of H in the general case: 

H= I EAPe 
This equation has the following meaning. For any ip Cz H and any function /, continuous on (t{H), 

where 

^i4E) ^ ^l^{{-cx,,E]) ^ WPEi^W^ 

This function increases monotonically with E and is upper semicontinuous (in jumps takes on the 
higher value). It fully determines the spectral measure /i^. In the above equation the integration 
is done with this measure. 

Here we open a parenthesis on measures. The main line of the discussion continues in Section 

5.8. 

5.6 Measures 

Let /i(x) be a real monotonically increasing upper semicontinuous function on R, finite at any 
finite x. With this function is associated a measure on the Borcl sets A of R: 



/z(A) := / d^(x) = / XA{x)d^{a 



which is nothing else than the total variation (increase) of /i on A, written as the Stieltjes-Lebesgue 
integral of XA(a;) with respect to nix). (The notation should not confound the reader: fi with 
a real number argument means the function, with a set argument the measure. In this way, 
/x(x) = /i((— oo,a;]) and fi{{x}) is the measure of the one-point set {x}.) The set function thus 
obtained satisfies, indeed, the properties of (positive) measures, namely, it is nonnegativc, vanishing 
on the empty set and countably additive: 

oo 

MU-iA,) = ^MA.) 

if the Ai are pairwise disjoint sets. 

Conversely, a positive measure fi which is finite on the semi-infinite intervals (— oo,x] can be 
used to define a monotonically increasing upper semicontinuous function, denoted also by /x, by 
setting 

fi{x) :=^((-oo,x]) 

This is called the distribution function of the measure fi. 

One can define the derivative of the measure fi with respect to the Lebesgue measure at the 
point X as 

- — (x) — lim/i(J)/TO(J) < oo 
dm Jlx 

provided the limit exists. The limit is taken on open intervals containing x and shrinking to x. 
If the distribution function is continuous and differentiable at x (perhaps with infinite derivative) 
then its derivative n'ix) = {dfi/dm){x). Oppositely, if {dfi/d'm){x) finitely exists then fJ.{x) is 
continuous and differentiable at x and n'{x) — {dfj,/dm){x). In the points of discontinuity of fi{x), 
fj,'{x) does not exist, but {dfj,/dm){x) exists and is infinite. The set 

S = {x : {dji/dm){x) does not exist finitely or infinitely} 

is a part of the set where n'{x) does not exist. Therefore we can apply a theorem ([124], Ch. IV, 
Theorem 9.1), valid for functions which are locally of bounded variation (like /i(x)), to conclude 
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that both the /i and the Lebesgue (m-) measure of S is vanishing, so (d/i/drn) (a;) is well-defined 
apart from a set of zero fi- and Lebesgue measure. 

Let ^ be a Borel set. We say that a measure /i is concentrated on A if /i(A^) — (A'^ = R\ A, 
the complement of A). The smallest closed set on which /i is concentrated is the support of fi, 
supp /i. So if ^ is concentrated on A then supp /i C A, the closure of A. 

An a; G R is a point of increase of the distribution function /i if for every £ > 0, fi{x + e) > 
fi{x - e). 

Problem 1. supp/i is the set of points of increase of ^{x). 

A Borel set A is called an essential support of jjl if /i(A'^) = and for any Aq <Z A such that 
/u(^o) = 0; 'm{Ai^)) = as well. So the support is unique and closed, the essential support is unique 
only apart from sets of zero ji- and Lebesgue measure and is normally not closed. Countable unions 
and intersections of essential supports are essential supports. An essential support may not be a 
subset of the support. An example is Ayj B where A is an essential support and B a nonempty 
Borel set with m{B) = and _Bn supp/i = 0. However, if A is an essential support then ^n supp/i 
is again an essential support, the closure of which is the support. 

It follows from the definitions, that any measure /i has an essential support on which (d/i/dm) (x) 
exists (maybe infinite). 

Any measure /i can be decomposed into pure point and continuous part, /i = /ipp + /icont- These 
are respectively characterized by the equations 

A^PpW = "^ l^ppiixt}) 

i 

that is, /ipp is concentrated on a countable set of points, and 

Mcont({a;}) = lim/i([a;-e,.T + e]) = 0, 

the /icont-measure of every a; G R is zero. If the function fi is continuous, the corresponding 
measure is purely continuous; if it makes jumps of heights a^ in the points Xi, 

d/ipp(a;) = 2_, d-i^ix — Xi)dx (5-6) 

i 

For /ipp there exists a smallest essential support: the set of all points of discontinuity of the function 
/i. 

If jjL and v are two measures, we say that /i is absolutely continuous with respect to v, and 
write /i <C i^, if /^(A) = whenever t^(A) = 0. In words, if v is concentrated on a set, then /i is 
concentrated on the same set. If jjl and v vanish on the same Borel sets, i.e., they are mutually 
absolutely continuous, we call them equivalent, and write ^ ^ v. Equivalent measures have the 
same support and essential supports. Most often v = m, the Lebesgue measure. If we say only 
that /i is absolutely continuous, we mean it with respect to the Lebesgue measure. Clearly, if /i is 
absolutely continuous, it is continuous. 

Problem 2. If two measures are absolutely continuous and have a common essential support 
then they are equivalent. 

Problem 3. Find an absolutely continuous and a pure point measure whose supports arc not 
essential supports. Find a measure, whose support is the smallest essential support. 

A function / defined on R is said to be absolutely continuous on a finite or infinite interval J if 
to every e > there exists a (5 > such that for any finite collection of disjoint open subintervals 
(ai,6i), (02,62), •••,(«»,&«) of J 

^\f{bi)- f{ai)\<e whcne-vcY ^{b,-ai)<5 (5.7) 

i i 

Since 5 does not depend on n, the property holds also if the number of intervals is infinite. 

Problem, 4- A measure is absolutely continuous if and only if its distribution function is abso- 
lutely continuous on every finite interval. 
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= {xeB 


ac 


= {xeB 


sing 


= {xeB 


sc 


= {xeB 


PP 


= {xeB 



We say that two measures /i and v are mutually singular, and write /i _L i^, if there are disjoint 
Borel sets A and B such that /i is concentrated on A and v is concentrated on B. Most often 
V — m and we simply say that ^ is singular. Any pure point measure is obviously singular. 

Problem 5. Any essential support of a singular measure is of zero Lebesgue measure. 

A measure is said to be singular continuous (with respect to the Lebesgue measure) if it is 
singular and continuous. An example is given in the next subsection. Thus, /x has the following 
decompositions into pairwise mutually singular measures: 

/^ /^PP ^ Mcont /^ac I /^sing /^ac ^" Msc i /^pp 

Problem 6. Show that any measure uniquely determines its pp, ac and sc parts. 

According to the value of {A^/drn) (x) , one can find essential supports for the different compo- 
nents as follows (taken over from Gilbert and Pearson [52]). Let 

i? = {a; e R : {dji / dm){x) exists, < (d/i/dm) (x) < oo} 

The sets M, Mac, Afging, Afgc, Afpp defined below are essential supports respectively for /i, /iac, Msing, /^sc, Mpp: 

< (d^/dm)(x) < oo} 

< (dy^/dTO)(x) = ^'{x) < oo} 

{d^/dm)(x) = oo} 

(d/i/d?Ti) (x) = oo, /x({x}) = 0} 

(d^/dm)(x) = oo, fJ,i{x}) > 0} 

Obviously, Mac,Msc and Mpp are pairwise disjoint sets. They are subsets of supp/i and Mpp C 
supp/ipp (Mpp is the smallest essential support of /Xpp), but Mac, Mging and Mgc may not be 
parts of the supports of the respective measures. As mentioned earlier, by intersecting them with 
the corresponding supports, we still get essential supports. It may happen, however, that Mac is 
nonempty but /Xac = 0: surely this is the case if ?TT,(Mac) = 0; or Mgc is nonempty, even uncountable, 
but fisc = 0. For more details see [52]. 

5.7 Cantor function 

The Cantor function is a continuous monotonically increasing function which grows from to 1 
exclusively in the points of the middle-thirds Cantor set 

oo 

C = {^ 3-"x„ : x„ = or 2} 

n=l 

In the interval [0,1] it is given by the formula 

oo oo 

a(x) = ^2-"-ix„ if x = ^3-"x„ , x„e{0,l,2} 

n— 1 n— 1 

with the remark that triadic rationals are taken with their infinite representation. More interesting, 
it is entirely determined by the conditions 

1. 2a{x) = a(3x) , < x < 1/3 

2. a{x) +a(l -x) = 1 

3. a{x) is monotonically increasing 

4. a{x) = if X < and a{x) = 1 if x > 1 

from which the moments or the Fourier transform of a can be computed. For example, the latter 
is 



e**^da(x)=e**/2j| cost/3" 

n=l 
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The function is represented by a 'devil's staircase'. The measure a is singular continuous, because 
it is continuous and its support, C, is of zero Lebesgue measure. Observe that the Fourier transform 
does not tend to zero when \t\ goes to infinity (check it with t — 27r3'^) as it would do if a were 
absolutely continuous. 

The support of a singular continuous measure may have positive Lebesgue measure. Define, 
for example, an odd function {p(—x) = —^{x)) on R by setting 

oo oo 

^i(x)^^2-'i^2-Pa{x-p/q) if a; > 

q=l j,=0 

This is a continuous function increasing in every point of R: /i(a; + e) — ^^'{x) > for every real x 
and positive e. ^{x) is therefore the distribution function of a continuous measure jjl whose support 
is R. Furthermore, UreQ(C + r) is a Borel set (a countable union of closed sets) of zero Lebesgue 
measure and is an essential support for ^; thus, ^ is singular continuous. This situation may occur 
for Schrodinger operators with unbounded potentials. If V^ = tanna and a is a Liouville number 
(an irrational number which is extremely well approximated by rationals), all the spectral measures 
are purely singular continuous, and there are spectral measures the support of which is R ([130]). 
Problem 1. Give an example of a pure point and a singular continuous measure having a 
common essential support. 

5.8 Spectral measures and spectral types 

In the comparison of two measures, /i and v, yielding /i <C i^, /i ^ J^ or /i _L i/, either or both can 
be projection- valued, and a projection- valued measure has support, essential support and uniquely 
determined pp, ac and sc parts, like real measures. In particular, for any vector ip, /i^ <C P, where 
P is the spectral resolution of the identity. Similar relations for the pp, ac, sc parts will be written 
down below. Let us start with the decomposition of the spectral resolution of the identity: 

^ pp ' -'cont -'ac I -'sing PP ' -'ac i -fsc 

Is it true that the different terms are spectral projection- valued measures? The answer is yes, but 
for the moment it is not even clear that Pac(A) and Psc(A) project onto _ff-invariant subspaces. 
By definition, 

Ppp(A)^ ^ Pi{E})=P{Ana,,iH)) 
BeAncrpp(-ff) 

which is a spectral projection. The smallest essential support of Ppp is app{H). Notice that 

Ppp(A) - P{<jpp{H))P{A) = [HppMA) 
where the two projections commute. Let ip G "Hpp. Then 

fi^{K) =. uf =. iiPpp(R)^ir = E ^'A(m) 

i.e., /i^ is a pure point measure and /i^ ^ Ppp. 
Also, by definition, for any Borel set A, 

Peo„t(A) = P(A)- Yl P{{E})=P{A\app(H)) 
BeAn<Tpp(H) 

is a spectral projection. The support of Pcont is crcont(-ff). Notice that 

Pcont(A) = (/ - [npp])P{A) = [Hcont]P{A) = P(app (F) = )P( A) 

where the two projections commute. Let ip G "Hcont- Then for any P G R 

^.^{{E}) ^ ||P({P})^f = ||P({P})[Hco„t]V'f - ||Pcont({P})^|P = 
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i.e., /i^ is purely continuous and /x^ ^ fcont- 

A purely continuous spectral measure can still be decomposed into absolutely and singular 
continuous parts. We would like to perform the analogous decomposition on "Hcont, crcont(-ff) and 

-fcont • 

Problems. 

1. ^^ is concentrated on a Borel set A if and only if P{A)ip — ip. As a consequence, li ip G 
P{A)n, supp/x^ C A 

2. HHip <C /x^. In particular, if /i^ is absolutely or singular continuous then fXHii is absolutely 
or singular continuous, respectively. (Remark. We know already that if /i^ is pure point, 
fiHi: is also pure point because xp € Hpp and HHpp C "Hpp.) 

3. If /i^ = (/i^)pp + (/ii/;)ac + (Mi/')sc then there exist orthogonal vectors f/'PP, ■0'*'^ and ip^^ such 
that = ip'P'P + ip^'^ + ip^"^ and (Af^)pp = /i^pp, {Ht(,)a.c — IJ-ip"" and (/x^)sc = /i^s^. ffint. Choose 
disjoint essential supports and apply 1. 

4. Let V* be a sequence of vectors converging to a vector ^p (i.e., ||-0* — tp\\ — > 0). Suppose that 
fi^i are singular measures. Then /z^ is singular, ffint. (i) ^^i -^ ^^ on Borel sets, (ii) Let 
Ai be an essential support to ^^i^ then UA^ is an essential support to fi^. 

It follows from Problem 2 that the vectors generating a purely ac spectral measure form an H- 
invariant subspace "Hac! similarly, the vectors giving rise to a purely sc spectral measure form an 
TJ- invariant subspace "Hsc- (That the vectors of the same pure spectral type form a subspace is 
the subject of Problem 7 below.) According to Problem 3, Hac and Hsc arc orthogonal and span 
'Hcont- Thus "Hac and "Hsc are necessarily closed and 

iLcont = /lac ffi rise 

(Problem 4 shows also explicitly that Hsc is closed.) Now we can consider the restrictions of H to 
"Hac and to "Hsc: these are selfadjoint operators, their spectra are a^dH) and <7sc{H), respectively. 
Absolutely continuous and singular continuous spectral measures are concentrated respectively on 
<Ts.c{H) and (Tsc{H). 

In summary, we obtained the following decompositions: 



Correspondingly, for any ■!/; G "H, 
and 



iL — r^pp ® iLa.c ffi iLs 



iP^tP'Pi' + tp''" + tp^- 



What about the decomposition of Pcont? Because 

["Hcont] ~ ["Hac] + ["Hsc] 

where the two projections are orthogonal, we get 

Pcont(A) == [Hcont]P(A) - [Hac]P(A) + [Hsc]P(A) 

with commuting operators in the products. 

Problem 5. Show that the projection-valued measures \Ha,c\P and ["Hsc]^ arc absolutely and 
singular continuous, respectively. 

Since the decomposition of Pcont into ac and sc parts is unique, we obtained that 

Pac=[Hac]P , Psc-[Hsc]P (5.8) 

So Pac(A) and Psc(A) project onto iJ- invariant subspaces within "Hac and "Hsc, respectively. It 
remains to show that they are spectral projections. 
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Let A and B be disjoint essential supports for Pac and Pgc, respectively, which arc disjoint also 
from (Tpp(_ff). Then 

[n^c] ^ P.c{A) = P{A) and [nsc] = Psc{B) = P{B) 

Indeed, for any ip £ H, ['Hac]V' = V'^^ while 

p{A)i, = p{A)r'' + p{A){r' + ^pp) - Pac(A)^ + (PsM) + Ppp(A))V' 

The proof is similar for [Hsc]- The result shows that Pac (A) and Psc(A) are the spectral projections 
P{A n A) and P{B n A), respectively. Using Eq.(5.8) and the definition of ^^, it is immediately 
seen that fj,^ <€. Pac if V" S "^ac and ^^ <C Psc if V' € "^sc- 
At last, 

a{H) ^ app{H) U aaciH) U (J,e(i?) 

where 



(5.9) 



It is also important to know, how to reconstruct the spectrum from real measures. Clearly, if 
/i is any real measure equivalent to P then in Eq.(5.9) P can be replaced by /i. 
Problem 6. Let /x be a measure and n^ <^ ii for all tp gH. Then P <C yU- 
Problem 7. fia-^i+i3^2 ^ fi^i + fi^2. 
Problem, 8. Let ■0^,-0^,... be a normalized basis in % and q > for z = 1,2,... such that 



a(7J) = 


= supp P 


'^pp(^) = 


= ess. supp Ppp 


<J.c{H) - 


= supp Pac 


<Tsc{H) - 


= supp Psc 



Z^i^i Ci < oo. Then 



P-H 



il^ip^ 



i=l 



In the particular case oi H — P{Ti) and H the Schrodinger operator, we have a much simpler 
result. Let 5^ and 5^ be the unit vectors concentrated on and 1, respectively. Then 

P- ^50+^^1 (5.10) 

Call /i the measure on the right side, /x <C P is obvious, one has to show P <C /i. Let A be such 
that ^(A) = 0. Then P(A)(5" = P{/\)5'^ = and, as a consequence, 

P{A)H"5° = P(A)i?"(5i =0 for n = 1, 2, ... 

From this we can conclude that P(A) ~ because of the following. 

Problem 9. The set of vectors {H'^d°,H"d^}'^^g is a basis in ^^(Z). 
Clearly, Eq.(5.10) holds if and 1 are replaced by any two successive integers. 

The sum of spectral measures is, in general, not a spectral measure. 

Problem 10. /i^i^^2 = /i^i + /i^2 if and only if the two vectors are in orthogonal i/-invariant 
subspaces. 

Therefore, for the Schrodinger operator the right member of the relation (5.10) is not a spectral 
measure, in general (but it is a spectral measure for the free Laplacian, see Section 5.9). 

If H has a pure point spectrum, and for the basis in Problem 8 we choose an orthonormal set 
of eigenvectors, ^Ciji^i is equal to the spectral measure belonging to ^aiV"*, where \ai\'^ — Ci. 
It is interesting to remark, that there always exist spectral measures equivalent to P, even if the 
spectrum is not pure point. The construction is suggested by Problem 10. Let A be an i/-invariant 
subspace. A vector ip^A'isa cyclic vector in A if any ^p £ A can be obtained as 



N 

im 

N 



lim VaAr„(i7"V-) 



n=l 
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In this case, ji^p <C /x^ for all (/? G .4, as one can see from Problems 2 and 7. Not all i/-invariant 
subspaces contain cyclic vectors (a counterexample is a subspace belonging to a degenerate eigen- 
value), but % can always be written as 

n = ®A , (5.11) 

a finite or infinite direct sum of orthogonal iJ-invariant subspaces with cyclic vectors. Let i/i' be a 
normalized cyclic vector in the zth subspace and cti 7^ 0, ^ jctip < oo. Any Lp eH can be written 
as (/3 = J2 f^iV^ where </?* <E Ai- According to Problem 10, 

^J■v = Msftv'' = Yl l/^»I^M>p' < Yl \^i\'^H' =" Msc.^' 

and, according to Problem 6, P ^ MsoiV'*' 

The decomposition (5.11) is nonunique. We can start with any tp^, choose Ai as the smallest 
subspace containing H^ip^ for every nonnegative integer n, choose any i/i^ orthogonal to Ai, and 
so on... . The smallest number of subspaces we need for the decomposition may be called the 
multiplicity of the spectrum of H . Problem 9 shows that the multiplicity of the spectrum of a 
Schrodinger operator on £^(Z) is 1 or 2. If the multiplicity is 1, that is, there exists a cyclic vector 
(in H), H is multiplicity free also in the sense defined in 5.5. 

Multiplicity may not be uniform on the spectrum, therefore it is appropriate to define it also 
locally: Let A be a Borel set such that P(A) 7^ 0, and Ha = P{A)H. The decomposition (5.11) can 
be performed on Ha , and the smallest number of subspaces we need to it is called the multiplicity 
of the spectrum on A. If A is an eigenvalue, we obtain the usual notion of multiplicity. 

5.9 A spectral measure for Hq 

Let H = e^{Z), H = Ho = T + T-^ (recaU: T is the left shift). We are going to compute exphcitly 
ligo . The way to proceed is to find a measure /i such that 

((5°, i/"(5°) = f E'^d^iE) for aU n 

Then we can identify /i with ^^0 . 
For any unitary operator U 

Choose for U the Fourier transformation which maps £'^(Z) onto L^^^ [0, 1]: For any ip G £'^(Z) 

00 
(C/V')(x)- J2 ^„e2-"- , a;e[0,l] 

n— — 00 

In particular, {US^){x) = 1 and UT^^U^^ correspond respectively to the multiplication by 
exp(=F27rzx); therefore UHU~^ corresponds to the multiplication by 2cos27ra:;. This gives 



1 - 2 

((5°,iJ"(5°)= f (2cos27ra;)"da; = 2 /"' (2 cos27rx)"dx = f E" 

Jo Jo J-2 



dE 



1-2 7rV4 - E^ 
where we applied the substitution E = 2 cos 2ttx. From here we recognize 



, ,„. dE/nV^^'W iiEe{-2,2) 

[ otherwise 

By integration, 

fj.so{E) = - + - arcsin(£;/2) if -2 < £■ < 2 

2 TT 

and below —2 and 1 above 2. This function is manifestly absolutely continuous, thus H^ has an 
absolutely continuous spectrum in [—2,2]. Because of the shift-invariance of Hq, ^gi — ^go and 
thus P ^ Id so- This shows that cr(iJo) = [^2, 2] and the spectrum is purely absolutely continuous. 
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Notice that (5° is not a cyclic vector: H"S'^, n — 1, 2, ... generate only the 'even' subspace, 

£ ~ {ip : tl)-k — i^k for all fc} 

This subspace is i^o-invariant and the orthogonal projection onto it is not a spectral projection 
(cannot be written as a function of _ffo only). Therefore _ffo is not multiplicity free. In fact, the 
spectrum is uniformly of multiplicity 2. 

5.10 £2(Z) versus ^^(N) 

Consider the matrix (5.1) of the Schrodinger operator on ^^(Z) and replace the matrix elements 
H^ifl = ifo,-i = 1 by zeros. What kind of effect this tiny modification can have on the spectrum 
of the modified operator H'l According to the'classical' Weyl theorem ([115], Vol. IV), 

acss{H') = cross (-ff) 

that is, apart from isolated eigenvalues, the two spectra coincide. The matrix of H' is of block- 
diagonal form, so H' can be written as H' = Hi (B Hr, where Hi and Hr are selfadjoint operators 
acting on ^^(-N \ {0}) and ^^(N), respectively. Thus, 

a{H')=<j{Hi)Ua{Hr) 

For strictly ergodic potentials acss{Hi) — acss{Hr), and, hence, = acss{H). 

The fate of spectral types under such a 'finite-rank' perturbation is partly uncertain. Point and 
singular continuous spectra may change into each other (Gordon [55], del Rio, Makarov, Simon 
[122]). However, the absolutely continuous spectrum is robust, 

(T.c{H) = a,,{Hi) U (Tac(ff.) (5.12) 

and this holds, whatever be the (selfadjoint) boundary condition at 0, defining Hi and H^- The 
reason is that the singular spectrum is sensitive to the boundary condition at 0, while the absolutely 
continuous spectrum can be fully characterized by the large-time behaviour of propagating wave 
packets and, hence, is insensitive to local perturbations. In fact, it is stable under more general 
(trace class) perturbations (see e.g. [77], Chapter X, Theorem 4.4). 

5.11 Asymptotic behaviour of generalized eigenfunctions. Subordinacy 

We started Section 5 with a remark about the sufficiency, for Physics, of polynomially bounded 
solutions. This remark is based on the following deep theorem: 

The spectral resolution of the identity has an essential support £ such that for every E ^ £ there 
exists a solution ip of the Schrodinger equation satisfying the following condition: For any 6 > d/A 
(d is the space dimension), with a suitably chosen positive constant c = c{S,E), 

\'ip{x)\ < c{l + \x\^)^ for ah a; in R'' or Z'' (5.13) 

The most general form of this theorem (going beyond Schrodinger operators) can be found in 
Berezanskii's book [22]. The Schrodinger case in the continuum, with precise conditions, assertions 
and proof is described in Simon's survey paper [128] (also in [29]). A simple proof for the i?^(Z) 
case, due to Lacroix, is presented in the book of Bougerol and Lacroix [25]. 

Clearly, the same result is valid separately for -Ppp, -Pac and Psc and is nontrivial for Pac and 
Psc- It, however, does not permit to distinguish between these two cases. 

It is expected that knowing the whole family of solutions for a given energy E, one can decide 
whether E is in the spectrum and which spectral type it belongs to. The right notion to deal with 
this question is subordinacy. It was introduced by Gilbert and Pearson [52] for the Schrodinger 
differential equation on the half-line, extended by Gilbert [51] to the problem on R, and by Khan 
and Pearson [78] to the discrete equation on N. A discussion of it can be found in Pearson's book 
[108]. Results for the equation on Z can be safely deduced from these works. 
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For a two-sided complex-valued sequence ip £^nd an integer N let 

N 



l^lk = {J2* \^n 



|2n1/2 



The star means that for A^ < the summation goes over n = 0, —1, ..., N. Fix E. A solution tp of 
the Schrodinger equation with energy E is said to be subordinate at +oo (— oo) if for every linearly 
independent solution (p, 

UWn/Mn^o 

as iV — > +00 (— oo). Clearly, there can be at most one subordinate solution at either side, and it 
suffices to check the condition on a single ip linearly independent of ip. Subordinate solutions for 
real energies are real. 

Below, by 'solution' we mean a solution of the Schrodinger equation for given E. Consider the 
following sets. 

• M' = R \ {E G R : there exist two solutions, one which is subordinate at +oo but not 
at — cx), the other which is subordinate at — oo but not at +oo}. 

• M'^^ — {E e R : every solution is nonsubordinate at +00} 

U{E (z R : every solution is nonsubordinate at —00} 

• Mg'; = {E e R : there exists a solution which is subordinate at ±00} 

• Mg(, = {E G R : there exists a solution which is subordinate at ±00 but is not in 

• M' = {E e R : there exists a solution which is subordinate at ±00 and is in i?^(Z)} 

Then M' is an essential support for P and M' = (Tpp{H). {M' is the disjoint union of M!^^, M^^ 
and M' and is a subset of a{H).) Furthermore, Pac 7^ if and only if m{M'^^) > 0, in which 
case M'^^ is an essential support for Pac (and so is M'^^ n CTac(-ff))- If Psc 7^ 0, M'^^ is an essential 
support for Pgc (and so is M'^^ n asc{H)). If J is a real interval, Pac{J) = Ppp{J) = and M!.^ n J 
is an uncountable set, then Psc (J) > 0. 

The above characterization shows that physicists are not quite right when saying that in the 
singular continuous spectrum 'all the solutions are critical'. In fact, there is one solution (the 
subordinate one) which is less (or more) 'critical' than the others and shows a clear analogy with 
an eigenvector. Notice, however, that a subordinate solution may not decay to zero at ±00. 

The union giving M'^^ corresponds to the union (5.12). On the intersection of the two parts of 
M!^^ the multiplicity of the spectrum is 2, otherwise it is 1 (Kac, [75], [76]). 

Recently, Al-Naggar and Pearson [3] developed further the characterization of the absolutely 
continuous spectrum in the case of the differential equation on the half-line. Due to Eq.(5.12), 
their result applies to the €^(Z) case in the form presented below. 

For fixed real E a complex solution -0 of the Schrodinger equation with energy E is called 
rotating at -|-oo (—00) if 

N N 

n=0 n=0 

as iV — > +00 (—00). Notice that a complex solution for real E is the linear combination of two 
real solutions with complex coefficients. Let 

M"j. ~ {E e R : there is a solution which is rotating at least at one of ±00}. 

Then, if m{M!l^) > 0, Pac 7^ and M"^, is inside an essential support of Pac- 

Unfortunately, it is not known so far whether M"^, itself is an essential support for Pac- 
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6 Schrodinger equation with strictly ergodic potentials 

6.1 Strict ergodicity 

The first step is to check minimality and unique ergodicity of the sequence defining the potential. 
We discuss this question only briefly; the reader may consult some other lectures of this School 
([1], [40], [95], [114]) and references therein. The minimality of a sequence (and then that of the 
hull) is equivalent to the infinite repetition, with arbitrary precision and bounded gaps, of every 
finite segment. In the case of sequences taking values from a finite set, once minimality is known, 
unique ergodicity means that the word frequencies exist (the defining limits converge). 

The cosine potential of the Almost-Mathieu equation is a uniformly almost periodic sequence: 
For any e there exists a sequence of integers, {rii}, with bounded gaps such that for all i and all n 

|K+n, -Vn\<e 

(see e.g. [24]). This implies also strict ergodicity. A weaker form of almost periodicity (in the mean- 
or Besicovitch sense, see [24] , which is equivalent to the existence of an atomic Fourier transform of 
the sequence, with square summable coefficients) also implies strict ergodicity. Sturmian sequences, 
or more generally, sequences generated by the circle map with A being the union of a finite number 
of intervals of the type [a, 6) (cf. point (2) in Sec. 2) belong to this class. The Fourier coefficients of 
the 1-periodic function Xa{x) are easy to compute. It is worth noticing that in this case the Fourier 
series converges everywhere but may not represent the sequence in a finite number of points. If 
this occurs, it is the function Xa, and not the series, which defines a strictly ergodic sequence. 

Some substitutional sequences, which are not Sturmian, also admit an atomic Fourier transform. 
Examples are the regular paper-folding and the period-doubling sequences. As one-sided sequences, 
these are strictly ergodic. 

In general, for every primitive substitution ^ one can build up two-sided minimal sequences, 
which are also uniquely (and, hence, strictly) ergodic ([56], [102]; see also [39]). The construction 
goes by concatenating a left- and a right-sided fixed point of some power of ^. More precisely, one 
can find two letters a and h and an n > 1 such that ^"(a) — ...a, C"(&) — b... and, with rj = ^", 
both u — ri°°{a) = ...a and v = 77°° (6) =6... contain the word ab and the two-sided sequence uv 
is minimal and uniquely ergodic. Among others, this holds for the Thuc-Morse and the Rudin- 
Shapiro sequences, although their Fourier transform is not atomic. Different constructions of a 
two-sided sequence (for example, via symmetric extension) may violate minimality (the starting 
sequence may not be recurrent). If both sides are related to the same substitution, the essential 
spectrum will not suffer, but the singular spectral measures can be seriously perturbed. (Compare 
with Section 5.10.) 

6.2 The spectrum of H{u) = Hq + V{u) 

Let us start with three remarks. 

1. The shift T is defined on £'^{Z) by {Tip)n — ipn+i- This is a unitary operator, therefore 
T"_ffr~" is unitary equivalent to and, thus, has the same spectrum as H for any finite n. This 
obviously holds for any potential. Notice that THoT^^ = Hq and TV{uj)T-'^ = V{Tu}), thus 

T''H{lu)T-'' = HiT'uj) 

2. If V'' is a sequence of potentials converging pointwise to a bounded potential V, then 
H'' = Ho + V" tends strongly to H = Hq + V, i.e., for any fixed V e ^^(Z), H'^'iJj -^ Hip. Indeed, 

||(ij'= - H)^j;\\' = Y. \^n - Kl'IV'nl' ^0 as fc ^ 00 

n 

3. If the bounded selfadjoint operators H'' strongly converge to the bounded selfadjoint H 
and A is an open interval not intersecting a{H^) for any fc, the spectral projections xa{H'') also 
strongly converge to xa{H) ([77], Theorem VIII. 1.15); on the other hand, they all vanish and, 
hence, xa{H) = 0. This, however, means that A does not intersect <j{H). In short. 



[ja{H'^)Da{H) 
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Problem 1. If V{ijj) is minimal, a{H{uj)) is independent of ui. Hint. Use the above remarks to 
prove that for any w and w', <j{H{lij)) D a{H{u!')) and a{H{uj')) D a{H{u!)). 

The cj-independent spectrum is denoted by (j{H). 

Minimality also implies that there is no isolated point in the spectrum. Even less is sufficient: 

Problem 2. If the potential is bounded and recurrent then (7{H) — crcss(-ff)- Hint. If E is an 
eigenvalue and tp the corresponding eigenvector, ip" = T^^^ifj with suitably chosen almost-periods 
fc„ is a Weyl sequence. 

We emphasize that the proofs of the above problems do not use ergodicity, only minimality or 
recurrence. This is the more interesting, because the same statements were verified for p-a.e. w by 
using ergodicity without minimality, see Section 5.2 for the random case. 

For one-dimensional ergodic potentials Pastur proved that any E is p-almost surely not an 
eigenvalue (Sec. 5.2). In fact, for this holding true, we need less than ergodicity: If fi is a shift- 
invariant set of potentials and p is a shift-invariant probability measure on J7, which is defined on 
cylinder sets, any fixed E is not in app{H [uj)) with p-probability 1. Indeed, 

Tx{E}{H{uj))T-^ =X{e}{H{Tlo)) 

and 



p oo ^ 

Y,J{S'',XiE}(H(T--Lo))S^)dp{co)= Yl y(<5°,Xm(HM)<5°)dpM 



n— — oo 



All the integrals exist, the first integral takes value in [0,1] {x{e}{H{u!)) = or it projects to a 
one-dimensional subspace of £^(Z)) and the last sum gives or oo. The value of the first integral 
is therefore 0, implying X{E}iH{u!)) = for p-a.e. w. 

The content of the above statement is that app(H(u!)) changes when oj changes. On the other 
hand, the closure of app{H{uj)) may not change: 
For ergodic potentials there are closed sets Spp, Sac and Sgc such that 



(Tpp{H{u})) = Spp 

asc{H{uj)) = Ssc 

for p-a.e. oj (Kunz, Souillard [91]). This theorem is the analogue of the assertion of Problem 1, but 
is a great deal more subtle than the constancy of the spectrum. Moreover, the result surely does 
not hold for all w, even if the potential is strictly ergodic: a counterexample of [71] was mentioned 
in Section 5.2. 

6.3 Integrated density of states 

For finite systems the integrated density of states (IDS) as a function of E counts the number 
of eigenvalues per unit volume below E. Let 7J^ £ denote the restriction of H{(jj) to the interval 
[—L, L] with Dirichlet boundary condition. The IDS for H^ £ is 

Nu),l{E) = X ( number of eigenvalues < E oi H^ £) (6.2) 

For uniquely ergodic potentials Mu:^l has an w-independent limit, Af{E), when L goes to infinity. 
As a limit of monotonically increasing functions, Af{E) is monotonically increasing and can be 
shown to be continuous. It is, therefore, the distribution function of a continuous measure on R, 
denoted also by J\f. Furthermore, it can be shown that suppA/" = cr{H). The use of Dirichlet 
boundary condition in the construction is not exclusive: any other boundary condition yielding a 
Hermitian restriction of H{uj) leads to the same IDS. (This, however, may not be true in higher 
dimensional spaces.) Obviously, the eigenvalues have to be counted with multiplicity. 
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The IDS for uniquely ergodic potentials can also be obtained by dealing directly with the 
infinite system. Fix any ^(a;) and let ^isn{E), n € Z, be the distribution functions of the spectral 
measures for H{lu), associated with the canonical basis in £'^{Z). Their average 

n— — L 

which is the distribution function of a measure {Ml), converges to M{E) as L goes to infinity ([8]). 

A part of the above results is easy to understand. If E is outside the spectrum, /o(w) = /x^o (E) 
is a continuous function on the hull and, because of unique ergodicity, its average along trajectories 
exists and is independent of w: this is J\f{E). If the spectrum is a Cantor set, the spectral gaps are 
dense everywhere (see Section 5.3) and M{E) has a unique extension from the gaps to R into an 
upper semicontinuous function. This, of course, does not explain, why is the spectrum a Cantor set 
and why is J\f{E) continuous. The IDS is constant in the gaps with values taken from a well-defined 
set (see Section 6.9). 

Recall that the (pointwise) convergence of the distribution functions does not imply the con- 
vergence of the measures on every Borel set. For example, if H{uj) has a pure point spectrum, for 
each E e crpp(i7(a;)), 

Ml{{E}) < 1/(2L + 1) ^ = MilE}) 

but A/'L(crpp(-ff(a;))) == 1 for each L while Af{<Tpp{H{Lu))) = 0. In general, the infinite summation, 
which occurs because of the cr-additivity oi Ml, does not commute with the limit i — > oo. 

The above remark makes less surprising the observation, that the type of the measure Af may 
have nothing to do with the type of the spectral measures. For example, in the case of random 
potentials, for p-a.e. uj the spectrum is pure point while H{lu) generates the same, continuous, IDS. 
In particular, Af is absolutely continuous for smooth p. For Bernoulli distribution and sufficiently 
large potential strength Af contains a singular continuous component with a support of positive 
Lebesgue measure ([28]). If the spectrum is a Cantor set of zero Lebesgue measure, a continuous 
Af is necessarily singular continuous: examples are presented in Section 7. If the IDS is singular 
continuous, the (differential) density of states can be given no meaningful definition. 

6.4 IDS and Lyapunov exponent 

In Section 3.2 we introduced the Lyapunov exponent, more precisely, jj.{E, V) and ^ ±{E, V) for 
an arbitrary potential V. According to a theorem by Fiirstenberg and Kesten [47], in the case of 
ergodic potentials there exists a function j{E) such that for every fixed E, for p-a.e. w the four 
numbers coincide to give 

7{E,io)= hm -^ln\\Ti^^a{E,io)\\ 

|n|^oo \n\ 

and "f{E,Lu) = 7{E). A detailed study of many related questions, as, for example, the uniformicity 
in E of the convergence can be found in Goldsheid's work [53]. 
Let 

S = {{E,Lo) : 7(i<^,a;) does not exist or 7^ 7(i<^)} 
Se = {lo : {E,Lo) e S} 
S'^ = {E:{E,uj)eS} (6.3) 

According to [47], for every E, p{Se) ~ 0. By the Fubini theorem, {m x p){S) = and, for p-a.e. 
w, m{S'^) = 0. 

In the theorem of Fiirstenberg and Kesten the restriction to p-a.e. uj is essential, even if the 
potential is strictly ergodic. A counterexample is provided by the Almost-Mathieu equation for 
A > 2 and a a Liouville number, where for any E G cr(i?) the set Se, defined above, is nonempty 
([8]). However, for potentials generated by primitive substitutions, Hof [63] proved that 5* is empty. 

There is a remarkable formula connecting the Lyapunov exponent to the IDS, found by Herbert 
and Jones [58] and Thouless [137]: 

-i{E) = I \n\E - E'\m{E') (6.4) 
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The main observation leading to this formula is the following. For fixed w let H^-^ L denote the 
restriction of H{uj) to the interval [1, L] with Dirichlet boundary condition. It has an L x L matrix 
with V{uj)i, ..., V{lo)l in the diagonal, 1 everywhere above and below the diagonal and elsewhere. 
Let '0 be the solution of H{uj)'>p ~ Eip with initial condition ipQ — 0, -ipi — 1. Then 

dct{E - Hlf^l^) = ^/jL+i 

Indeed, both members are polynomials of E of degree L, they have the same roots and the same 
principal coefficient. Let Ei^^, ...,El^l be the roots, then 

L 

L-i In \iJL+i{E)\ = L-i Y.^n\E- E,^l\ 

i=l 

When L -^ oo, the two sides go to the respective sides of Eq.(6.4), apart from a set of {E,uj) 
of zero m X p measure. Equation (6.4) holds, in fact, for all E £ C. The proof is based on the 
subharmonicity of "/(E) and is due to Herman [59] and Craig and Simon [35]. 

6.5 Results on the set {E : -f{E) = 0} 

Let us recall that for periodic potentials 

(TiH) - a,,{H) = Go where Go ^ {E : 7(S) = 0} 

For Ho, the pure kinetic energy operator, we still have a{Ho) = [—2, 2]. 
For almost periodic potentials Deift and Simon [38] found 

m(Go) < 4 

with equality if and only if the potential is constant. In fact, this result extends to any ergodic 
potential (Last [92]). 

Ishii [68] and Pastur [106] obtained the following result: 
Let V{uj) be a bounded ergodic potential and P" denote the spectral resolution of the identity for 
H{uj). Suppose that m(Go) > 0. Then for p-a.e. w, P^ is concentrated on Gq. 
The proof is easy: Choose a typical uj, such that m{S'^) = (cf. Eq.(6.3)). The complement of 
Go can be written as Gq = AU B, where 

A = {E:-f{E,uj)^j{E)>0} , B = Gg n 5"^ 

Now P^(i?) — because P^ is absolutely continuous and B is of zero Lebesgue measure. On 
the other hand, for every E Cz A, either there is no polynomially bounded solution or there is an 
exponentially localized solution. Therefore P!^^{A) = and, hence, Pi;^(Go) = 0. 

There is a different way to formulate this result. We can drop out of Go any set of zero Lebesgue 
measure and close the rest: P^^ is still concentrated on this set which is now closed and, therefore, 
contains <Tac{H{uj)). A particular closed set to which this remark applies is the essential closure of 
Go. 

The essential closure of a set A is 

3°"" = {E : m{A n (E - e, E + e)) > for all e > 0} 

A short inspection may convince the reader that A'^^^ is indeed a closed set, it is inside the closure 
of A, and what it does not contain from A is a set of zero Lebesgue measure: 

m{A\A'''')^0 

Thus, the Ishii-Pastur theorem can be brought into the form 

a,,{H{uj)) CG^ 
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holding for p-a.e. lu. 

What do we drop from Gq when we take its essential closure? Imagine that, apart from an 
absolutely continuous spectrum, H has also a singular continuous spectrum on a disjoint set D of 
zero Lebesgue measure. If 7(i<^) = on _D, D belongs to Go but not to its essential closure. Similar 
is valid for a point spectrum on which the eigenvectors are not exponentially localized. 

Most remarkably, the converse of the Ishii-Pastur theorem is also true and we have the following. 
For p-a.e. uj 

aac(i/M)=Gr (6.5) 

Moreover, Go is an essential support for P^^. 

This theorem was proven by Kotani for the differential equation ([86], [87]) and adapted by 
Simon to the difference equation ([129]). 

Remarks. 
(i) o-a.c{H{uj)) is the support of Pl^^, see Eq. (5.9), but it may not be an essential support. In fact, 
what we add to Go when taking its essential closure, may be a set of positive Lebesgue measure, 
see the example of Problem 5.6.3. 
(ii) Since (j,^c{H{u))) n Go is also an essential support for Pi;|^, we have, in particular, 

m(crac(i:f(a;))nGo) = m(Go) 

for p-a.e. lo. 

(iii) The main content of the above theorems is that the absolutely continuous spectrum of H{uj) 
is p-almost surely nonempty if and only if ?7i(Go) > 0. The 'if is clear from the preceding remark; 
the 'only if holds because m(Go) — implies that the essential closure of Go is empty. 

6.6 The role of periodic approximants 

Given a bounded aperiodic potential V, one can define a sequence of periodic approximants, V'^ , 
in such a way that V^ converges to V pointwise. According to the second remark made in Section 
6.2, H^ = Hq + V'^ then tends strongly to H — Hq + V and Eq.(6.1) holds for the spectra. Our 
aim is to minimize the covering set on the left side of this equation, by choosing the best periodic 
approximants. 

The best periodic approximants of potentials of the type V{uj)n — g{na+ui)^ where 5 is a period- 
1 function and a is irrational, are obtained by replacing a with its best rational approximants, Uk- 
For example, if a = (v5 — l)/2, ak = Fk-i/Fk, where Fk is the fcth Fibonacci number. 

If V is a substitutional potential, V'' can be chosen to be the periodic repetition of the sequence 
f{oJn), evaluated on the word rf{a)rf'{b) (cf Sections 2 and 6.1). 

6.7 Gordon-type theorems 

Minimality sometimes implies that the potential repeats itself (exactly or with very good precision) 
on three neighboring intervals, one of which starting with 1, and this holds for an increasing 
sequence of interval lengths. In such cases, it can be shown that the Schrodinger equation has no 
solution decaying at the infinity, and thus the spectrum is purely continuous. The first proof of 
this kind was given by Gordon [54]. Suppose that 

(y_L+i...y_iVb) = {Vi...Vl-iVl) = {Vl+i...V2L-iV2l) (6.6) 

Then 

and, applying the two sides of the Caley-Hamilton equation 

AI-(ItAl)Al + I = ^ (6.7) 

to the vectors ^_l and ^Pq (cf Eq.(3.1)), one finds that for any solution t\} of the Schrodinger 
equation 

max{||*_L||,||*L||,||*2L||}> ^ll^oll 
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Clearly, if (6.6) holds for an increasing sequence L„, no solution can decay on both sides, so there 

can be no eigenvector to H . Let us see two examples. 

(i) Avron, Simon [8]: For the Almost-Mathicu equation, for A > 2 and a Liouville number (for 

every positive integer k there are integers pk and qt such that \a — Pk/qk\ < fc '"')j for a.e. uj the 

spectrum is purely continuous. In this case, Eq.(6.6) does not hold exactly but with an extremely 

large precision. The spectrum is, in fact, purely singular continuous: this comes from the Ishii- 

Pastur-Kotani theorem, because j{E) > for all E ([4]). 

(ii) Delyon, Pctritis [43]: For circle potentials (see Section 2) with A being an interval, for every 

A, a.e. Lo and a.e. a the spectrum is purely continuous. In order that Eq.(6.6) hold true, some 

weak condition on the continued fraction expansion of a has to be imposed. This limits the result 

to almost every a. 

A variant of the Gordon theorem uses only two intervals. Suppose that in Eq.(6.6) the second 
equality is verified for an increasing sequence L„. This situation arises, for example, with substi- 
tutional potentials, if the substitutional sequence starts with a square. Applying the two sides of 
Eq.(6.7) to ^0, one obtains 

max{|trAL|||*L||,||*2L||}> ^ll^oll 
This yields the absence of decaying solutions at +oo for energies such that 

is a bounded sequence. In some cases one can show that this set of energies is just the spectrum 
(see [135] and Section 7). 

6.8 Kotani theorem for potentials of finite range 

Any periodic potential on £'^{Z) is of finite range, i.e., takes on a finite number of different values. 
Any Schrodinger operator on £^(Z) with a periodic potential has a purely absolutely continuous 
spectrum. Therefore, the following theorem by Kotani [88] may be surprising. Let V{uj) be an 
ergodic nonperiodic finite-range potential, p the corresponding measure on the hull of V . Then, 
for p-a.c. Lu, H{u)) has a purely singular spectrum. 

This theorem is at the origin of many results on Schrodinger operators with circle- and substitu- 
tional potentials; we discuss them in the next section. It has long been an open question, whether 
the restriction to p-a.e. lo can be dropped, and whether the spectrum is purely singular continuous 
for every u, if the potential is strictly ergodic. A recent result by Hof, Knill and Simon [64] goes 
in this direction: The authors show that for strictly ergodic potentials, either H{lu) has pure point 
spectrum for all lo or there is an uncountable dense set, although of zero p-measure, in the hull 
for which the spectrum is purely singular continuous. This latter case is shown to be realized for 
circle potentials if a is irrational and A is a half-open interval. This includes Sturmian potentials. 
Furthermore, for potentials generated by primitive substitutions, the Lyapunov exponent ^{E,lo) 
is independent of w (Hof, [63]); Kotani's theorem then implies that the spectrum is purely singular 
for all LO. 

6.9 Gap labeling 

Gap labeling is a book-keeping for spectral gaps. The IDS naturally assigns a number to each 
gap: the (constant) value Af{E) for E in the gap. The question is then to characterize this set of 
numbers. The first example of gap labeling for an almost periodic Schrodinger equation was given 
by Johnson and Moser [73] in the continuum case. The method used in the discrete case is quite 
different. 

Problem 1. Let V be an L-periodic potential and suppose that H = Hq + V has no missing 
gap. Show that N{E) — k/ L, k — 1,2, ..., L — 1 in the gaps. 

For strictly ergodic potentials, the values of the IDS in the gaps are taken from (but do not 
necessarily exhaust) a set S, determined by an algebraic theory, the K-theory. The application of 
K-theory to gap labeling was developed by Bellissard, Lima, Testard, Bovier and Ghez ([19], [16]). 
A detailed discussion can be found in [13], [14]. 
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Sometimes it is possible to describe S quite explicitly. For example, according to Bellissard, 
Bovier and Ghez ([14], [13], [16], [27]), for circle and substitutional potentials S is the module 
(smallest additive group) containing the word frequencies of V{uj), restricted to the interval [0, 1). 
This module contains the integers (the sum of the frequencies of cylinders with common base (n 
and k in Eq.(2.4)) is 1), therefore S* is a group with respect to addition modulo 1, like for periodic 
potentials. 

For Sturmian potentials this gives 

S = {ka + m(l - a) : fc, m e Z} n [0, 1) = {ka + to : fc, m G Z} n [0, 1) 
Notice that for a — K/L, where K and L are relatively primes, this gives 

12 L-1 



in accordance with the assertion of the Problem above and the example (ii) in Section 4. If a is 
irrational, S is countably infinite and dense in [0, 1]. If there exist, indeed, gaps corresponding to 
values in a subset S" of S which is still dense in [0, 1], the spectrum is necessarily a Cantor set and 
the IDS is continuous. 

In our discussion of the periodic Schrodinger equation (Section 4) we introduced the notion of 
a missing gap (or closure of a gap). When there are no intervals in the spectrum, as is often the 
case with strictly ergodic potentials, what one can unambiguously assert is the absence of missing 
gaps or the completeness of gap labeling. Surely, there is no missing gap if 

RaugA/" = {M{E) : i<^ is in a spectral gap} — S , (6.8) 

the set of all the admissible values. (Notice that each admissible value is taken on in at most a 
single gap.) As in the periodic case, there is no general method to check whether or not all the 
gaps open. This is known to hold true in a few cases, like the period doubling potential or the 
Fibonacci potential with A > 4. 

Example. For the classical Cantor set, in the /cth gap on the nth level (n = 0, 1, ...; fc = 1, ..., 2") 
the Cantor function a{x) — {2k — l)/2"+^. These numbers form an additive group modulo 1. The 
same values are taken by the IDS of the hierarchical Hamiltonian ( [90] , see also [89] ) in the domain 
of the parameters where the potential is limit periodic. In this case the gap labeling is complete. 

7 Schrodinger equation with Sturmian and substitutional 
potentials 

7.1 Fibonacci potential 

Fibonacci substitution was the first to be used to define a two-valued potential (see Section 2) and 
to study the spectral problem of the corresponding Schrodinger operator (Kohmoto, Kadanoff and 
Tang, Ostlund et al. [79], [80], [81], Casdagh [30], Siito [135], [136]). This is also the first example 
where the trace map ([2], [84], [111]) was fully exploited. 

The Fibonacci sequence, as any minimal sequence, has almost-periods: these are the Fibonacci 
numbers F„ (Fo — Fi — l,F„+i = Fn + i^,i-i). Choose, for example, 

Vn = Xilin + l)a] - [na\) (7.1) 

then the almost-periodicity is expressed by the equations 

Vi+Fr, = T^ if n > 3 and 1 < ? < i^„ 
V^„j.,„ = V^ if n > 1 and 1 < ? < i^2«+i 

The transfer matrices over the almost-periods replace the single transfer matrix Ti^l in the L- 
periodic case. For an L-periodic potential 

Ti^2L = T^^L (7-2) 
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Let 

For the Fibonacci potential Eq.(7.2) is replaced by 

This, together with dct M„ = 1, implies for r„ = tr M„ the recurrence relation 

Tn+2 = Tn+lTn — T„_i 

The initial conditions (with the choice (7.1)) 

T_i = 2 , tq=E , Ti= E- X 

contain only 2 parameters, and this indicates that the recurrence must have a nontrivial invariant. 
Indeed, 

2,2,2 A r n2 \2 

"^n+l +'''«+ Tn-1 ~ Tn+lTnTn-i - 4 = (Ti - TqJ = A 

From the trace map and the invariant, the following results can be deduced ([135], [136]). 

1. a(H) = {E : {t„(^)}^^i is bounded} 

(Remember that for an L-periodic potential <t{H) = {E : \trTi^L{E)\ < 2}.) 

2. The spectrum oi H(uj — 0) is purely continuous (the proof is of Gordon- type with two 
intervals) . 

3. j{E, u; = 0) = for aU E e cr(H). 

4. The Lebesgue measure m{(j{H)) ~ (found by confronting the Ishii-Pastur-Kotani theorem 
with Kotani's theorem for potentials of finite range). As a consequence, iy{H) is a Cantor 
set (because (y{H) contains no isolated point). 

5. The spectrum of H{uj = 0) is purely singular continuous. 

6. For all E e (j{H) all solutions of H{ix: — O)'0 — Eip are polynomially bounded (lochum, 
Raymond, Testard [66], [67]). 

It is an open question whether the spectrum is purely singular continuous for all uj. Let us recall 
that TH{uj)T^^ — H{Toj) = H{uj + a) is unitary equivalent to H{lu), therefore the singular 
continuity holds true for a; = ka (mod 1), k € Z, which is a countable dense set in [0, 1]. Due to 
[64], this result has recently been extended to an uncountable dense set of uj, still of zero measure 
(cf. Section 6.8). 

7.2 General Sturmian potentials 

The case of Sturmian potentials for arbitrary irrational a was investigated by Bellissard, lochum, 
Scoppola and Testard [17]. The methods applied to the Fibonacci case can be extended to treat 
the general problem, even though this is technically more involved. All the results found for the 
Fibonacci potential remain valid and the same questions are unanswered. 

There is an interesting numerical work by Ostlund and Kim [OK] and a rigorous study by 
Bellissard, lochum and Testard [18] on the a-dependence of the spectrum oi H = Hq + V^"\ 
where 

cf. Section 2. For rational a the potential is periodic and the spectrum is the union of a finite 
number of intervals. When approaching a rational a from above and from below, the numerical 
plot of the energies belonging to the spectrum clearly reveals a discontinuity, see Figure 2. [18] 
show that this reflects the discontinuity of the characteristic function which generates the potential. 
The gap edges vary continuously on irrational a's. For a rational value r, 

liml/(") 9^1imy(") 
35 



and none of them equals y^'"-'. In fact, the two hmits are not periodic, only ultimately periodic 
(as if we took a periodic sequence, cut off a finite segment which is not a period and glue the two 
infinite pieces together). Both yield the same essential spectrum as Ho + V^'^'' and both create 
(different) isolated eigenvalues in the gaps. 

7.3 Period doubling potential 

The potential is generated by the period doubling substitution ^(a) — ab, ^(6) ~ aa. The corre- 
sponding Schrodinger equation was studied by Bellissard, Bovier and Ghez [15], [27]. The trace 
map is a fundamental tool of the analysis. Results 1.-5. valid for the Fibonacci potential hold true: 
The spectrum is a Cantor set of zero Lebesgue measure, it is purely singular continuous (again, not 
known for all elements of the hull) and the Lyapunov exponent vanishes in the spectrum. There is 
a rather detailed knowledge about the gaps, their behavior as a function of the potential strength. 
Gap labeling is complete: all admissible values are taken by the IDS in the gaps. 

7.4 Thue-Morse potential 

Historically, the first studies were done by Axel et al. on the phonon frequency spectrum of the 
harmonic chain with masses generated by the Thue-Morse substitution ^(a) = ab, ^(b) = ba. 
In [9], [10], [11] we find the proof that the phonon spectrum is a Cantor set, a numerical work 
on its box dimension suggesting zero Lebesgue measure, gap labeling and a study of generalized 
eigenfunctions. This latter established the existence of extended states for a dense set in the 
phonon spectrum. 

The Schrodinger equation with the potential generated by this substitution has also been widely 
studied, namely by Bellissard [12], Delyon and Peyriere [43] and Bovier and Ghez [26], [27]. Valu- 
able numerical work was done by Riklund et al. [121]. Bellissard [12] identified the set RaUgA/" 
(cf. (6.8)) and studied the dependence of the gap widths on the potential strength. Delyon and 
Peyriere [43] proved the absence of decaying solutions and, hence, the continuity of the spectrum. 
They showed that the generalized eigenfunctions are not 'too small' on a geometric progression. 
The continuity of the spectrum was obtained as a byproduct also in [15], [27]. Bovier and Ghez 
[26], [27] proved, in a more general context (see next section), that the spectrum is of zero Lebesgue 
measure. With the continuity this implies that the spectrum is purely singular continuous. 

7.5 Systematic study of substitutional potentials 

Bovier and Ghez [26], [27] succeeded to find a general condition on primitive substitutions, assuring 
that the corresponding Schrodinger operator has a spectrum of zero Lebesgue measure. Their work 
deeply exploits the Kolar-Nori results on trace maps ([84], see also [111])- 

The action of a substitution ^ can be defined on transfer matrices and on traces of transfer 
matrices. ^ acting on traces is called a trace map. ^ carries a letter a of the alphabet A into a word, 
the transfer matrix corresponding to a into a product, i.e. a monomial, of transfer matrices in 
reversed order, and the trace of the transfer matrix into a polynomial of traces of transfer matrices 
belonging to letters of an enlarged alphabet B. These generalized letters are special words over 
the original alphabet A, not containing repeated letters of A. Therefore B is also finite, and one 
can see that the trace map is closed in the following sense: 

^ carries the trace belonging to any letter of B into a polynomial of traces belonging to letters of 
B. 

Let bi, ...,bq be the letters of B and Xi a real variable associated with bi. The idea of Bovier 
and Ghez is to retain the highest degree monomial 

Q 
f^{xl,...Xg) = Yl^j' 

i=i 

of the polynomial corresponding to bi and to define a substitution cj) on B such that 4'{bi) contains 
bj kij times. Based on earlier experience, namely with the Fibonacci substitution, one may expect 
that, by imposing some conditions on (j), it is possible to control the high iterates of the trace map 
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and, hopefully, the spectrum. The right notion is semi-primitivity, a property which is somewhat 
weaker than primitivity (see [26] for the definition) and is straightforward to verify. The authors 
prove that, if ^ is primitive, (f) is semi-primitive and the substitutional sequence contains the word 
bb for some b (z B, <j{H) is a set of zero Lebesgue measure. Although the proof is more complicated 
than in the cases discussed earlier, the main idea is again to show that the Lyapunov exponent 
vanishes in the spectrum (the existence of the word bb is used here) and, then, to confront the Ishii- 
Pastur-Kotani theorem with the Kotani theorem for finite-ranged potentials. If the substitutional 
sequence starts with a square, the two-interval version of the Gordon theorem yields also that the 
spectrum is purely singular continuous. 

The conditions of the theorem are fulfilled by many named substitutional sequences as, for 
example, the Fibonacci, Thue-Morse, period doubling, binary non-Pisot, ternary non-Pisot and 
circle sequences. Apart from the binary non-Pisot sequence, singular continuity is also verified in 
the cases listed above (but it is not proved for the other potentials in the hull) . A notable exception 
is the Rudin-Shapiro substitution which is primitive and therefore the spectrum is purely singular 
([64]), but for which the substitution (f) is not semi-primitive and, hence, it is not known whether 
lyiH) is of zero Lebesgue measure. 

8 Solutions to the problems 

2.1 

The proof is based on the observation that for bounded sequences pointwise convergence is 
equivalent to convergence in the metric d. 

(i) 1. implies 2. 

Indeed, s eVL(s) = VL(t). 

(ii) 2. implies 3. 

Q,{s) is bounded: for all t £ ^(s), d{s,t) < 6sup|s„|. 

Suppose that s is not almost-periodic. Then there exists an e > and an infinite sequence of 
real intervals {ak,bk) such that 6^ — a^ — > oo and d{T"'s,s) > e for every integer n G U{ak,bk). 
Choose an integer Uk in (afe, bk) such that Uk — ak-^oo and bk — Uk — >■ oo. Since fl(s) is closed and 
bounded (compact in the topology generated by d), T"''s has at least one limit point t in U{s). 
Fix an arbitrary integer n. We have the triangle inequality 

d{s,T''t) > d{s,T'"'+"s) - d{T'"'+"s,T"t) . 

For k sufficiently large, Uk + n E {ak,bk) and thus d{s,T"''^'^s) > e. On the other hand, 
d{T"''s,t) -^ implies that diT^^+^s^T^t) ^ as /c ^ oo. It follows that d{s,T"t) > e for 
all integers n. Hence, s cannot be in r2(t), which contradicts 2. 

(iii) 3. implies 2. 

Let t € r2(s), t = limj^oo T^^ s, where lim means pointwise limit. Fix an e > and let {nk} be 
the almost-periods with gaps smaller than £g. Given j, choose k such that nk-i < ij < rik- Then 
vij = nk — ij < is and 

d{s,T"'H) < d{s,T'''s) + d{T'''s,T"'H) 
< e + d{T"'^T'^s,T"^H) 

Take the limit j — > oo. Since < ttij < £g, there will be an jti G (0, i^) which occurs infinitely many 
times among the m^'s. Taking the limit only on the subsequence of j values for which rrij = m, 
we obtain d{s,T"^t) < e. Since e was arbitrary, this proves that s € n(t). 

(iv) 2. implies 1. 

Suppose t G i^{s). Then T"i e f7(s) for every integer n, because 0(s) is shift-invariant. J7(s) 
being closed, this implies fl{t) C U{s). On the other hand, according to 2., s € n{t). Interchanging 
s and t, the above argument yields f2(s) C n{t). 

2.2 

Let A C fl(s) be a T-invariant set. One has to show that p{A) ~ or 1. Suppose the opposite, 
i.e., < p{A) < 1. Define a probability measure /i by setting 

n{B)=p{AnB)/p{A) 
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It is easy to verify that /i is a T- invariant probability; on the other hand, it differs from p (e.g., on 
^{s) \ A). This contradicts the uniqueness of p. 

2.3 

Let w e fi. If {r"w}^Q is an infinite set, the sequence contains no repetition. Therefore 



1 > p({r"c.}- o) - E Ki^"^}) = E ^({^}) 



and hence p({w}) — 0. 
3.2.1 

The proof goes by induction. trA" = tr/ = 2. Muhiplying the Caley-Hamihon equation 

A^ - (trA)^ + / = 
with A^'"^ and taking the trace one finds 

ir A^ = iv Ail A''-^ - ii A^-'^ 

in which one may recognize the recurrence relation for the Chebyshev polynomials. Setting iv A 
2 cos a, the solution of the recurrence is 

iiA^ = 2 cos ka 



5.1.1 

Let Hijj — Eip where V' is subexponential: for any a > there exists a 6 > such that 

IV'fcl < &e"l'=l for aU k, 

but -0 is not in £'^{Z). It follows that either X^fcLo IV'fcP = oo or X]fc=-oo iV'fcP = oo. Suppose, for 
instance, the first. Define ^" — '0/(X]r=o iV'iP)"'^''^ fo^ n > 0. Then Hip"- — Eijj"-, and (p", given by 

n _l i'l if < fc < n 
'^ 1^ otherwise 

is a normalized vector for all n. We show that a suitable subsequence of {<p"}^i is a Weyl 
sequence. It is clear that for all fc, (^^ ^ with increasing n, so (p" goes to zero weakly. On the 
other hand, 

{Hip'')k^Eip'^ if fcT^ -l,0,n,n+l 

This implies that 

if^ = if fc7^-l,0,n,n + l 

(cf. Eq.(5.3)). It is easy to check that 

r" — ?/;" f" — — ?/;" F^ — —ih^ F^ — 7/;" 

=■-1 — 'A'O =-0 ^ 'P-l ^n — 'Pn+1 ^n+1 — Vn 

Therefore 

where the vectors on the right correspond to the definition (3.1). The first term goes to zero as n 
increases, the second, in general, not. However, there exists a subsequence Uk such that ||*5^^|| — > 
with increasing k. Indeed, if the opposite were true, one could find a positive constant c such that 
||*"|P > cfor ahn > 0. Now 



c < ii*;:ip = ii*«ii VE i^-^i' ^ 2||vi/„ii V E II* 



.11^ 



fc=0 k=0 
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and therefore 

n-l n-2 

*;=0 k=0 



2 



which contradicts the supposed subexponential nature of ip. The Weyl sequence we were looking 
for is {^"''j'^^i- 

5.2.1 

In general, let E be any accumulation point of the eigenvalues of H, say, En — > -E as n — > oo. 
Let "0" be the corresponding orthonormal eigenvectors. Then ip" form a Weyl sequence. First, 
they go to zero weakly: for any vector (p ^ H, the Bessel inequality 

oo 

ii^f >Ei(^"'^)i' 

Tl=l 

implies that {ip", (p) —¥ with increasing n. Second, 

Hi;" = EnV = E^" + {Er, - E)i;" 

and II (E„ - E)yj'^\\ ^ |E„ ~ E\ ^ 0. 

5.2.2 

Choose a sequence n^ such that Vn^. — > -E as fc — > oo. If (5" is the unit vector concentrated on 
n, 

\\{E-V)-'S"''\\^l/\E-VnJ^^, 

so the inverse of i? — F is unbounded. By definition, this means that E is in the spectrum. 
One can also construct a vector ip such that the solution oi {E — V)(p = ip for (/?, given by 

is not square-summable. Write E — cos 27r0, then trivially 

|^-K| < 27rmin|na-0-p| 

pGZ 

On the other hand, by the Kronecker theorem (Theorem 440 in [57]), there exists an increasing 
sequence Uk of positive integers and a sequence pk of integers such that 

\nka- -pk\ < 3/nfe 

The example is obtained by choosing ipn = E — Vn for n ^ Uk and zero otherwise: WipW^ < 6tt^ but 
(pn^ — 1 for all k, so (p is not square-summable. 

5.6.1 

If X is not in supp fi then for sufficiently small e > 0, fi{x + e) — fJ-ix — e) = ii{{x — e,x + e)) =0 
and, hence, x is not a point of increase of the distribution function. If x G supp /i then for every 
e>0, 

^i{x + s) — /i(x — e) = ^{{x — e, a; + e]) > /i((x — £, x + e)) > , 

that is, X is a point of increase. 

5.6.2 

Let /i, V be ac measures, A a common essential support. Suppose that ^{B) ~ 0. v{Br\A'^) = 
because A is an essential support of v. m.{B n A) = 0, because _B n A C A, an essential support of 
/x, and ^{B n A) = 0. But v is ac, so v{B n A) = 0. We found v <^ ji. The opposite is obtained 
by interchanging ^ and v. 

5.6.3 

(i) Let C be a 'thick' Cantor set (i.e., m{C) > 0), 

CC I lOO ^1C I lOO I ikn J 
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where J„fe are the open intervals appearing in the construction of C (cf. Sec. 5. 3). Let ,fnk{x) be a 
continuous function which is strictly positive on Jnk and vanishes outside J„fe, and 



^ ^ / /„fc(a;)dx < oo 



n=Ofe=l 

Then ^^ ^j, /„fe is the density of an absolutely continuous measure /i. supp/i = R (because the 
closure of C^ is R), and an essential support of /i is C"^. So /i(C) = but m{C) > 0, therefore 
supp /i is not an essential support of /i. 

(ii) For < p < q, p and g integers which are relatively primes, let fi{{p/q}) = 2^^^*. 
supp/i — [0, 1] while the rational numbers in (0, 1) form an essential support. /i([0, 1] \ Q) = but 
m([0, 1] \ Q) = 1, so the support is not an essential support. 

(iii) Let fj, be concentrated on a discrete set (whose only accumulation points can be ±oo), for 
instance, /x(R) — X]n=i A'd'"'})- Then supp/i is the smallest essential support. 

5.6.4 

(i) If iJ,{x) is absolutely continuous (ac) on any finite interval then the measure n is ac. 
Indeed, let A be a Borel set, m{A) = and suppose first that A is covered by a finite open interval 
J. Fix £ > and let J > correspond to e (cf. Eq.(5.7)). According to the definition of a set of 
zero Lebesgue measure, one can find an open set O such that A C O and 171(0) < S. O is the union 
of disjoint open intervals (ai,6i), (02,62), •••, all in J, so that A C U{ai,bi) and ^(6i — Oi) < S. 
Thus, 

m(^) < Yl A^((««' ^0) = XI (^(^*) " '"('^*)) < ^ 

Since e is arbitrary, /i(j4) = 0. 

If A is not covered by a finite interval, it is covered by countably many finite intervals Jk- 
fi(A n Jfc) = for all fc, thus fi{A) = 0. 

(ii) If /i is an ac measure, i.e., 11 <^ m, then ij.{x) is an ac function on any finite interval. 
We show this by proving that, given a finite interval J, for any e > there exists a (5 > such 
that for Borel sets A C J, m{A) < S implies /i(A) < e. Suppose this does not hold true. Then 
there exists an e > such that for all n > 1 and for suitably chosen A„ C J, m(An) < 1/2" but 
^i{A„) > e. Let B„ = U^„+iAfe. Now m(B„) < 1/2" and /z(S„) > /i(^n+i) > s. Since B„ C J 
for all n and i?„ is a decreasing sequence, it has a limit B C J for which m{B) = and ii{B) > e. 
This contradicts fi <ti m. 

5.6.5 

Suppose fi is a singular measure. Then there is some set A of zero Lebesgue measure such 
that /i is concentrated on A. The intersection of any essential support with A yields an essential 
support of zero Lebesgue measure. If B is an essential support and m{B) > then iJ.(B \ A) = 
and m{B \ A) > 0, a contradiction. 

5.6.6 

The definition of /ipp is constructive, see Eq.(5.6), therefore the pp part is unique. Let 

A^ = A*pp + Msc + Mac ^ A^PP + Msc + A*ac 

If /igj. 7^ /igj., there exists a set A such that to(A) ~ and nl^i^) 7^ /ig^(A). But then /iac(^) ¥= 
Mac(A) which contradicts /ii^CA) = Mac(A) = 0. 

5.7.1 

Let C be the middle-thirds Cantor set, a the sc measure whose distribution function is the 
Cantor function. Let fi he a pp measure whose smallest essential support is the set of boundary 
points {cni} of the complement of C, cf. Sec. 5.3. Then supp/i = supp a = C and C is also an 
essential support for both measures. 

5.8.1 

The following statements are equivalent: /i^ is concentrated on A, /i^(R) = /i^(A), UfAlP = 
||P(A)^||2, ||(/ - P{A))ij\\ = 0, (/ - P(A))^ = 0, P{A)iP = V- 

If ^|J £ P{A)H then P{A)tl; = ip, so fi^iA") = and supp/i^ C A. 

5.8.2 
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Let A be any Borcl set such that /i^ is concentrated on A. Then from Problem 5.8.1, P{A)ip = 
ip. H commutes with P{A), therefore P{A)Hip — Hip. From Problem 5.8.1 it follows that nn-ip is 
concentrated on A. By definition, we get fXHip ^ fJ'ip- 

If fit/j is singular, it is concentrated on a set A such that m{A) = 0. Therefore fXH^p is also 
singular. If /x^ is continuous, for any E eH, fi^{{E}) = \\P{{E})-)p\\'^ = 0, therefore P{{E})->p = 0. 
Hence, 

^^H4,{{E}) = \\P{{E})Hi;f = \\HP{{E})ijf = , 

that is, fiHip is continuous. So fxn^) is sc if /x^ is sc. If fi^ is ac then /iHi/" is ac because m{A) = 
implies /i^(^) = 0, which implies HH-ipiA) = 0. 

5.8.3 

Let A, B, C be disjoint essential supports respectively for (/x^)pp, (/i^)ac, (/^i/')sc- MiA i^ concen- 
trated on A U i? U C. From Problem 5.8.1, 

ip = P{A UBU C)ij = P{A)ij + P{B)'iP + P{C)i; = ip^^ + ip'''' + ip""" 

Now 

P{A)P{B) = P{A)P{C) = P{B)P{C) = 

therefore ip^^ , ip'^'^ , ip^"^ are pairwise orthogonal, and for any Borel set A 

fi^iA) = \\p{A)i^r = \\p{A)p{A)i^r + \\p{A)p{B)i^r + \\p{A)p{c)i^r 

= fi^pp (A) + fi^^c (A) + ^^sc (A) 
Using the relations P{A)P{A) = P{A n A), etc., we find 

/^^pp(A) = /A0(An A) == (^^)pp(A) 

and so on. 
5.8.4 

(i) Let B be any Borel set, then 

fi^.{B) = \\P{B)r\\' - (V.,P(B)^) + (^*-^,P(B)^*) + (^,P(i?)(#-V)) 
= fi^{B)+0{\\r-^\\)^^i4,{B) 

(ii) ip^ are singular, so there exist sets Ai such that /^^i(A?) = and ■m{Ai) — 0. Now 
m{UAi) = 0, for alH /x^i((UAj)'^) </i^i(A°) = andby (i), /z^((Uylj)^) = 0. So ^^ is concentrated 
on a set of zero Lebesgue measure and, hence, is singular. If all the fi^i are pp or all are sc then 
/x^ is pp or sc, respectively, because 'Hpp and 'Hcont are closed subspaces. 

5.8.5 

(i) [H-aclP is ac: Let i? be a Borel set, ra{B) = 0. Take any ip G H, then 

\\[nac]P{B)lPf = j|P(B)[Hac]^|P = A^^-(B) = (MV.)ac(S) = 

therefore [na^]P{B) = 0. 

(ii) ["HscJ-P < ["Hcontl-P — ^cont, therefore ['Hsc]^ is a continuous measure. We have to prove 
that it is singular. Choose an orthonormal basis {ipi}- Let Ai be sets of zero Lebesgue measure, 
^^sc concentrated on Ai. Take A — UAi. m{A) = and we show that \hLsc\P is concentrated on 

A.' 

||p(A^)[Hsc]^.|| = ll^(^^)^f II = ||P(n,A^=)^f II = |in^(^,')^(A^)^f II -0 

The last equality holds because PiyA'i)'^f = 0, according to Problem 5.8.1. It follows that 
[-Hsc]P(A'=)V'» = for ah i, and therefore \R^^\P{A'') = 0. 

5.8.6 

If \i{B) = then ii^{B) = for ah -0 e "H. Accordingly, ||P(S)V'|| = for all -0 e "H which 
means that P(-B) = 0. 

5.8.7 
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Let B be a Borel set, n^i{B) +/x^2(B) = 0. Then \\P{B)^P^\\ = \\P{B)i:^\\ ^ and hence 
P{B)ip^ = P{B)ip'^ = 0. This imphes 

f^c^^+priB) - |a|2||P(B)^i||2 + |/3|2||p(s)^2|j2 

+ a*(3{i;\P{B)i;^) + a(3*{4,^,P{B)^b^)=0 

5.8.8 

According to Problem 5.8.6, it suffices to show that fi^ ^ ^Cifi^t for ah ^ ^ %. Writing 

ip — ^ttiV'*: one proves 

exactly as in Problem 5.8.7. 

5.8.9 

(i) The set is linearly independent: Taking any finite linear combination with not all coefficients 
vanishing and expanding the sum in the canonical basis {(5"}, at least one 5" appears with nonzero 
coefficient. 

(ii) The set generates the canonical basis {S"}^^_^- For n = one obtains S'^,d^ and by in- 
duction one finds that S^", (5^"+^, ..., S", 5"+^ can be expanded with the vectors {H^5^ ^ H^^^yk=Q- 

5.8.10 

Let -0^ G T-Li and -0^ e 'H2 where T-Li are subspaces of "H. 

/x^i+^2(i?) = \\P{B)(^^+^^)f 

= ^l^^ (B) + ^i^2 (B) + (P(B)0\ V') + ii^^PiB)^') 

and 

{P{B)i;\i;^) = {P{B)[ni]ij\ [%]^') - 

for every Borel set B if and only if the two subspaces are orthogonal and one of them, say "Hi, is 
7J-invariant (and, hence, ["Hi] commutes with P{B) for all B). Then, however, H2 C Hi which is 
iJ- invariant and contains Tp"^. 

6.2.1 

Choose a sequence {uk} such that V'^ = V{T"''uj) -^ V{lu') pointwisc. The corresponding 
operators H'' = Hq + V^ all have the same spectrum, a{H{uj)), and converge strongly to H{uj'). 
Equation (6.1) holds in the form a{H{u!)) D <t{H{uj')). The opposite is also true because V{uj') is 
also minimal. 

6.2.2 

Any point in the spectrum which is not an eigenvalue is in the essential spectrum. So we 
have to prove that all the eigenvalues are in the essential spectrum. Let E be an eigenvalue, 
Ip the corresponding normalized eigenvector. Because V is recurrent, there exists an increasing 
sequence {nk} such that ¥>" = T^i^VT-"" (Vj; = K+nJ tends to V pointwise. Let (p^ ^ T-"*^-0. 
This goes to zero weakly ((cp, (^'^) — > for any fixed vector ip), so it is a Weyl sequence if e'' = 
{Ho + V - E)lp^ -)■ (in norm!). 

||gfe||2 ^ ||ynfc£fc||2 ^ ||(^^ + V''- E)ijf = \\{V'' - V)i^f 

There exists a sequence of integers nik — >■ 00 and a sequence rjk -^ such that \Vjl; — Vn\'^ < ijk ioi 



n 



< ruk, and J2\n\>mk l^"l^ < "^k- Therefore 



ll(^'-^)^lP = ( E + E )(^„'-K0'IV'nP<r/fe(l + 4sup|K. 

\n\<mk |n[>mfc 

SO e*^ goes to zero, indeed. 
6.9.1 
According to Section 6.3, the IDS can be obtained as 

N'iE) = lim MnLiE) 



42 



where MnL{E) is 1/nL times the number of eigenvalues less than E of i/„L, the restriction of H 
to an interval of length nL. A convenient choice of boundary condition is 

V'(A: + nL) = e*"VW 

The matrix of HnL has Vi, ..., VnL in the diagonal, 1 above and below it and e~*" and e*" in the 
upper right and lower left corner, respectively. Finding the eigenvalues of HnL is equivalent to 
finding the E values for which the transfer matrix Ti^nhiE) = Tl\^j^{E) has eigenvalues e^*" or 
trace 2 cos a: this can be seen from 

(cf. Section 3.2) and the boundary condition. Choose a ~ 7r/2, then one has to look at the 
zeros of trT{\.^(£'). From Figure 1 and Problem 3.2.1 one learns that this function has exactly n 
zeros in each band of the spectrum of H; therefore the IDS in the gap just above the fcth band is 
kn/riL = k/L, independently of n, which proves the assertion. 
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